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The  following  conditions  regulate  the  use  of  computer  programs 
developed  by  the  Hydrologic  Engineering  Center  (HEC),  Corps  of 
Engineers,  Department  of  the  Army. 

1.  The  computer  programs  are  furnished  by  the  Government  and  are 
accepted  and  used  by  the  recipient  individual  or  group  entity 
with  the  express  understanding  that  the  United  States  Government 
makes  no  warranties,  expressed  or  implied,  concerning  the 
accuracy,  completeness,  reliability,  usability,  or  suitability 
for  any  particular  purpose  of  the  information  or  data  contained 
in  the  programs,  or  furnished  in  connection  therewith,  and  that 
the  United  states  Government  shall  be  under  no  liability 
whatsoever  to  any  individual  or  group  entity  by  reason  of  any  use 
made  thereof. 

2.  The  programs  belong  to  the  United  States  Government. 
Therefore,  the  recipient  agrees  not  to  assert  any  proprietary 
rights  thereto  nor  to  represent  the  programs  to  anyone  as  other 
than  Government  programs . 

3.  The  recipient  may  impose  fees  on  clients  only  for  ordinary 
charges  for  applying  and  modifying  these  programs. 

4.  Should  the  recipient  make  any  modifications  to  the 
program(s),  the  HEC  must  be  informed  as  to  the  nature  and  extent 
of  those  modifications.  Recipients  who  modify  HEC  computer 
programs  assume  all  responsibility  for  problems  arising  from,  or 
related  to,  those  modifications.  User  support  from  the  HEC  to 
third  party  recipients  will  only  be  provided  after  the  second 
party  demonstrates  that  program  difficulties  were  not  caused  by 
their  modifications. 

5.  This  "Condition.-  of  Use"  statement  shall  be  furnished  to  all 
third  parties  that  receive  copies  of  HEC  programs  from  the 
recipient.  Third  party  recipients  must  be  notified  that  they 
will  not  receive  routine  program  updates,  correction  notices,  and 
other  program  services  from  the  HEC  unless  they  obtain  the 
program(s)  directly  from  the  HEC. 

6.  All  documents  and  reports  conveying  information  obtai-ed  as  a 
result  of  the  use  of  the  program(s)  by  the  recipient,  or  others, 
will  acknowledge  the  Hydrologic  Engineering  Center,  Corps  of 
Engineers,  Department  of  the  Army,  as  the  origin  of  the 
program (s ) . 
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HYDROPOWER  ANALYSIS  PROGRAM 
USING  STREAMFLOW  DURATION  PROCEDURES 

THE  HYDROLOGIC  ENGINFF.RTNG  CENTER 

1.  ORIGIN  OF  THE  PROGRAM 

This  program  was  developed  by  the  Hydrologic  Engineering  Center, 

U.S.  Army  Corps  of  Engineers.  It  is  based  on  techniques  used  for  the 
Corps  National  Hydropower  Study  in  1979  (HEC,  1979c).  Up-to-date 
information  and  copies  of  the  program  can  be  obtained  from  the  Center 
upon  request. 

2.  CAPABILITIES  OF  THE  PROGRAM 

The  HYDUR  program  has  been  designed  to  perform  hydroelectic  power 
potential  analysis  based  on  flow-duration  principles.  The  basic  pur¬ 
pose  of  the  program  is  to  integrate  a  flow-duration  (transformed  to 
capacity-duration)  relationship  to  determine  energy.  A  wide  range  of 
options  based  on  the.  flow  duration  concept  are  available  for  analysis. 

The  options  include  such  important  items  as  considering  the  number  and 
performance  characteristics  of  generating  equipment  (efficiency,  operating 
range  and  overload  capacity);  availability  of  water  and  head;  and  tailwater 
effects  on  generation.  The  program  is  thus  specifically  designed  to  analyze 
run-of-river  type  hydropower  projects  -  that  is  projects  whose  storage  is  not 
specifically  operated  to  alter  the  manner  of  hydropower  generation. 

The  program  has  two  options  for  providing  hydrologic  data.  A  flow 
duration  curve  may  be  directly  input  or  streamflow  data  may  be  input 
and  the  program  derived  the  flow-duration  curve.  Options  are  available 
to  further  adjust  streamflow  through  such  mechanisms  as  ratios,  losses 
and  diversions. 

The  program  can  also  produce  preliminary  estimates  of  hydroelectric 
power  potential  that  would  commonly  require  conventional  period  of  record 
sequential  analysis  such  as  is  provided  by  the  computer  program  HEC-5 
(HEC  1982) .  These  options  were  carried  over  from  the  analysis  procedures 
used  for  the  national  Hydropower  Study  in  which  the  goal  was  to  perform 
simple  type  analysis  for  virtually  all  potential  pi  jects.  These  options 
should  be  used  with  due  regard  for  their  inherent  assumptions  thus  guarding 
against  the  temptation  to  use  the  simpler  HYDUR  analysis  when  in  fact  more 
refined  procedures  are  warranted.  Options  that  fall  into  this  group  are 
those  that  permit  analysis  of  storage  projects  (reservoirs  operated  to 
conserve  runoff  for  later  release),  provide  for  preliminary  estimates  of 
construction  costs,  and  provide  for  national  level  valuation  of  power 
output.  A  few  of  the  automated  cost  and  benefit  items  may  be  overridden 
with  specific  input  data  and  more  capability  to  specify  site  specific  costs/ 
etc.,  are  planned  as  future  improvements  to  HYDUR. 
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The  power  benefit  values  contained  in  the  program  are  1978  regional 
capacity  and  energy  values  developed  by  the  Federal  Energy  Regulatory 
Commission  (FERC  1978).  The  cost  estimating  procedures  were  developed 
by  the  Corps  North  Pacific  Division  and  were  intended  for  use  in  the 
Nationally  scoped  National  Hydropower  Study. 

The  program  also  has  a  utility  feature  that  permits  automatic  sizing 
of  project  features.  The  program  can  select  the  optimum  capacity  of  a 
plant  based  on  maximizing  or  minimizing  any  of  several  criteria.  The  user 
is  cautioned  to  be  wary  of  optimum  capacities  that  are  determined  when  using 
the  cost  and  generalized  benefit  functions  available  in  the  program. 

Results  should  be  verified  for  reasonableness  and  are  appropriately  con¬ 
sidered  to  be  an  indicator  of  possible  plant  capacity  that  should  be  refined 
by  further  study  based  on  site  and  project  specific  information. 

3.  COMPUTER  REQUIREMENTS 

The  program  available  for  distribution  is  written  in  FORTRAN  IV, 
using  140K  words  of  central  memory  on  a  CDC  7600  computer  or  approximately 
1400K  bytes.  The  program  uses  four  scratch  units.  Unit  1  requires  an  89- 
character  record  for  each  input  card  that  is  read.  Unit  2  requires  a 
record  for  each  data  card  image  read  from  the  alternative  file;  the  record 
length  is  dependent  on  the  data  format  requested  by  the  user.  Unit  3  (Unit 
8  on  Harris  version)  is  used  in  conjunction  with  the  alternative  file. 

Unit  4  is  used  as  a  dump  file  which  contains  superfluous  output.  In  ad¬ 
dition,  to  being  used  as  a  separate  program,  HYDUR  can  be  linked  as  a  sub¬ 
routine  to  any  other  program  as  long  as  the  proper  variables  are  passed 
through  the  call  statement.  The  Corps  of  Engineers  version  on  the  Lawrence 
Berkeley  Laboratory  (LBL)  computer  system  has  additional  capabilities  that 
allow  it  to  read  flow-duration  information  directly  from  GETUSGS  data  files 
(IIEC,  1979a). 

4.  METHODS  OF  COMPUTATION 

a.  Determination  of  a  Flow-Duration  Curve 

The  flow-duration  curve  can  be  calculated  from  monthly  streamflows, 
daily  streamflows,  or  it  can  be  provided  directly.  In  the  first  two  cases 
the  program  counts  the  number  of  flow  values  less  than  or  equal  to  each  of 
70  discharge  intervals  defined  by: 

Q.=  10(l_1)/1°for  i  =  1,70  (Eq.l) 

An  additional  interval  count  is  made  of  the  flow  between  zero  and  1  cfs. 
This  equation  was  adopted  because  it  eliminated  the  need  to  sort  the  stream- 
flow  data.  The  logarithmically  spaced  intervals  provide  a  range  of  flow 
intervals  that  give  adequate  definition  of  the  flow-duration  curve  for  all 
size  streams. 

The  number  of  occurrences  greater  than  or  equal  to  the  flow  in  each 
interval  are  determined  and  divided  by  the  total  number  of  occurrences; 
thus  yielding  ordinates  that  describe  the  percent  of  time  each  interval 
is  exceeded.  The  ability  to  develop  seasonal  streamf low-duration  curves 
and  perform  subsequent  seasonal  power  analysis  is  only  available  in  this  mode, 
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When  the  curve  is  provided  directly,  at  least  15  to  20  points  should  be 
provided  to  define  the  curvature  because  linear  interpolation  is  used  to 
estimate  values  between  the  ordinates. 

Once  the  flow-duration  curve  has  been  developed,  consideration  to  the 
type  of  project  and  its  operation  must  be  evaluated  to  determine  if  additional 
adjustments  to  the  curve  are  warranted.  Analysis  of  a  run-of-river  project 
requires  no  adjustment  to  the  flow-duration  curve  since  any  inflow  to  the 
project  must  be  either  used  for  power  production  or  spilled  downstream. 

A  storage  project,  however,  is  capable  of  accumulating  excess  inflow  for 
future  use  during  low-flow  periods,  thereby  converting  the  inflow-duration 
curve  into  a  flatter  or  less  peaked  outflow-duration  curve.  If  the  input 
flow  duration  curve  or  flow  record  was  derived  from  regulated  flows  down¬ 
stream  from  the  project,  no  further  adjustment  is  needed.  If  the  flows 
are  input,  storage  should  be  accounted  for.  Analysis  can  be  performed 
externally  by  hand  or  machine  and  values  input  as  above,  or  if  a  simplified 
accounting  of  the  effects  of  storage  is  acceptable,  the  program  has  an 
option  that  will  perform  the  adjustments.  The  technique  developed  to 
perform  this  storage  effect  adjustment  is  discussed  in  Exhibit  2  entitled, 
Adjustment  of  Flow-Duration  Curve  for  Storage  Effect.  For  now,  however,  it 
is  important  to  note  that  once  an  adjustment  is  performed,  if  at  all,  the 
remaining  procedure  applies. 

Figure  1  below  illustrates  the  basic  analysis  procedure.  End  points 
are  established  at  0  and  100%  (Maximum  and  minimum  recorded  flows)  and 
analysis  proceeds. 

Power  is  generated  when  the  flows  on  the  streamflow  duration  curve 
are  less  than  QSUB,  (Submergence  flow),  the  flows  on  the  curve  are  greater 
than  QMIN,  the  minimum  flow  necessary  to  drive  a  turbine  and  the  operating 
head  is  above  minimum  operating  head  limits. 


Figure!.  Streamflow  Duration  Curve 
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b.  Calculation  of  a  Capacity  Versus  Percent  of  Time  Exceeded  Curve 


The  flow  duration  curve  is  transformed  to  a  capacity  duration  curve 
as  the  next  step  leading  to  energy  calculators. 

The  portion  of  the  flow-duration  curve  between  points  P^and  P^  in 
Figure  1  is  broken  into  40  evenly  spaced  percent  of  time  exceeded  intervals. 
Forty-one  discharge  ordinates  are  calculated  to  define  the  40  intervals  by 
linearly  interpolating  between  either  the  computed  or  the  user-supplied 
streamflow  duration  coordinates.  The  capacity  duration  curve  is  then 
constructed  from  these  41  coordinates  as  follows: 

1)  The  Power  Equation: 

(CAP)i  =  0.084603  *  (HW-TW)  *  EEF  *  QA^  (Eq.  2) 

subject  to:  QA_^  _<  QDES  (turbine  check)  (Eq.  3) 

( CAP )  _<  CAPDES  *  0VLOAD  (generator  check)  (Eq.  4) 

The  available  discharge  to  produce  power  is  defined  by: 

QAi  =  (Qi  *  QFACT)  -  DIV  -  QL0SST  for 

QSUB  <  Q  <  QMIN  or  (Eq.  5) 

QA  =  0  for  Q  <  QMIN  or  Q  >  QSUB. 

where:  (CAP)^  =  capacity  in  KW  based  on  the  flow  QA^; 

Q^  =  discharge  (cfs)  from  the  flow-duration  curve; 

QA.  =  discharge  (cfs)  of  flow  that  is  available  to 
produce  power; 

QSUB  =  Q  corresponding  to  tailwater  submergence; 

QMIN  =  Q  necessary  to  drive  turbine; 

QDES  =  maximum  penstock  capacity  in  cfs; 

HW  =  headwater  elevation  in  feet; 

TW  =  tailwater  elevation  in  feet; 

EFF  =  efficiency  (E^  *  E  *  E^)  expressed  as  a 
decimal  fraction;  ®  J 

CAPDES  =  installed  capacity  of  the  plant  in  kilowatts; 

0VL0AD  =  overload  factor  of  the  plant,  expressed  as  a 

decimal  (i.e.,  1.15  for  a  15%  allowable  overload); 

QFACT  =  user  supplied  factor  expressed  as  a  decimal  used  to 

adjust  the  streamflow  discharges  (e.g.,  drainage  area  ratio); 

DIV  =  flow  diversion  in  cfs  above  the  powerhouse.  (Average 
evaporation  losses  can  also  be  included  in  DIV.) 
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DIV  is  subtracted  from  the  streamflow  values  on  the 
flow-duration  curve  before  performing  the  power  analysis; 

QL0SST  =  diversion  of  water  in  cfs  around  the  powerhouse  (fish 
ladder,  leakage,  etc.)  This  flow  is  not  used  for  the 
power  production  but  effects  the  headwater  and  tailwater 
elevations ; 

E  =  efficiency  of  the  turbine  expressed  as  a  decimal  fraction; 

E  =  head  losses  of  the  overall  power  configuration  expressed 
®  as  a  decimal  fraction; 

i  =  one  of  41  percent  of  time  exceeded  ordinates; 

0.084603  =  62.4/737.56,  where  62.4  is  the  weight  in  pounds  of  1  cubic 
foot  of  water  at  50  F  and  737.56  is  the  conversion  factor 
from  1  kilowatt  to  737.56  ft-lbs/sec. 


2)  Additional  Hydrologic  Parameters 

The  amount  of  spillage  QS_.  =  (Q_.*QFACT)  -  DIV  -  QL0SST  -  QA_^  (Eq.  6) 

subject  to  QS  >0 
i  — 


Whenever  QA^  is  less  than  or  equal  to  QDES  and  is  also  within  the  flow 
limits  defined  by  QSUB  and  QMIN,  then  QS^  will  be  zero. 

The  HW  term  can  be  specified  as  a  function  of  total  reservoir  releases  QR^, 
which  is  expressed  as: 


QR.  =  QA.  +  QS.  (Eq.  7) 

The  efficiency  EFF  can  be  specified  as  a  function  of  QA.. 

The  TW  term  can  be  specified  as  a  function  of  the  tailwater  flow  QT.. 

QT^  is  defined  in  one  of  two  ways  depending  upon  user  selection: 

QT±  =  QAj,/UAPD  +  QSi  or  QT±  =  QA^UAPF  (Eq.  8) 

Where  UAPF  represents  a  factor  between  0  and  1,0  which  is  used  to 
adjust  the  daily  mean  flow  available  QA.  to  a  value  which  represents 
the  actual  flow  passing  through  the  penstocks  during  times  of  power 
operation  only. 

The  spillage  QS.  is  added  to  the  tailwater  flow  whenever  this  flow 
will  effect  power  operations.  In  configurations  where  the  power¬ 
house  is  remotely  situated  from  the  inpoundment  and  spillway  structure, 
the  spill  effect  should  not  be  included  (See  Figure  2). 
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Figure  2.  Schematic  of  Spill  Effect 


3)  The  program  makes  a  final  adjustment  of  the  capacity  duration 
curve  to  guarantee  that  it  monotonically  decreases.  A  capacity  curve 
that  does  not  so  decrease  is  not  logical.  The  adjustment  is  a  simple 
reordering  of  values  in  decreasing  order. 

c .  Calculation  of  Average  Energy 

Average  energy  is  calculated  by  integrating  the  area  under  the 
capacity  verses  percent  of  time  exceeded  relationship  derived  in  the 
previous  steps.  Figure  3  is  a  schematic  of  the  curve. 


Figure  3.  Capacity  versus  Percent  of  Time  Exceeded  Curve 


P2  C 

AE  =  F(8.76)  *p  \(CAP).  dp 


F(8.76)  * 


1 2 


i=P 


1 


Si  *(CAP).  *  AP 
"3  1 


(Eq.  9) 
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where:  AE 

( CAP ) j 

(dp, A  P) 


average  energy  in  megawatt -hours  (MWH); 

capacity  in  KW  at  the  ith  percent  of  time  exceeded 
ordinate; 

percent  of  time  exceeded  interval; 


Pi  =  lower  limit  of  integration  which  is  equal  to  the 
percent  of  time  exceeded  ordinate  associated  with 
QSUB; 

p2  =  upper  limit  of  integration  which  is  equal  to  the 
percent  of  time  exceeded  ordinate  associated  with 
QMIN; 

Sj  =  Simpson's  integration  constant,  where: 

S.  =  1  for  i  =  1,  41  (first  and  last  ordinates); 

=  4  for  i  =  2,4,6  .  .  .40  (even  ordinates); 

=  2  for  i  =  3,5,7  .  .  .  39  (odd  ordinates); 

8.76  =  8760  hours/1000  KW,  a  constant  that  converts 

KW  to  the  MWH  generated  in  a  year; 

F  =  length  of  season  expressed  in  years. 

(0.083  <_  F-<  1.0)  For  annual  analyses,  the  value  of 
F  is  1.0.  For  seasonal  analyses,  the  values  of 
F  must  be  greater  than  or  equal  to  1/12  and  less  than  1.0. 


d.  Calculation  of  the  Plant  Factor 


The  plant  factor  is  defined  as  the  ratio  of  the  average  load 
on  the  plant  for  a  given  time  period  to  the  aggregate  rating  of  all 
the  generating  equipment  installed  in  the  plant. 

PF  =  AE/  (CAPDES  *  8.76  *F)  (Eq.  10) 


where: 


PF  =  Plant  Factor;  (for  a  given  time  period); 


AE  =  Average  energy  in  MWH; 

CAPDES  =  Installed  plant  capacity  in  KW.  8.76  is  a  conversion 
factor  based  on  (1,000  KW/MW)  /  (8760  hours/year); 

P  =  1.0  for  annual  power  analyses. 
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e.  Determination  of  Dependable  Capacity  and  Annual  Firm  Energy 


The  HYDUR  program  includes  routines  for  determing  approximate 
estimates  of  dependable  capacity  and  annual  firm  energy.  The  conventional 
approach  for  usch  estimates  requires  sequential  analysis  of  power  production 
coupled  with  the  system  load  relationship.  Approximate  methods  based  on 
flow-duration  principles  of  analysis  have  been  developed  and  included  in  HYDUR. 

IACWR  defines  dependable  capacity  as  the  "load-carrying  ability 
of  a  station  or  system  under  adverse  conditions  for  the  time  interval  and 
period  specified  when  related  to  the  characteristics  of  the  load  to  be  supplied." 
Annual  firm  energy  is  defined  (IACWR,  1965)  as:  "electric  energy  which  is 
intended  to  have  assured  availability  to  the  customer  to  meet  all  or  any 
agreed  upon  portion  of  his  load  requirements." 

The  flow-duration  curve  constructed  from  the  reservoir  inflows 
does  not  account  for  the  reduction  of  the  inflow  variability  that  occurs 
when  a  reservoir  is  constructed.  It  is  quite  possible  that  while  reservoir 
inflows  may  be  zero  for  several  months  during  the  critical  period,  carry¬ 
over  storage  may  maintain  dependable  discharges  from  the  project.  Therefore, 
the  minimum  power  releases  from  a  project  are  generally  considerably  higher 
than  the  minimum  reservoir  inflows.  The  magnitude  of  this  effect  will 
depend  upon  the  variability  of  the  inflow  regime,  the  storage  capacity  of 
the  reservoir,  the  release  schedule  of  the  project,  and  the  length  of  the 
drought  period.  Thus  the  flow  duration  curve  to  be  used  for  power  analysis 
must  account  for  storage  effects.  If  a  storage  adjusted  curve  is  not  in¬ 
put  directly,  an  approximately  estimate  of  storage  effects  can  be  considered 
by  activating  the  HYDUR  storage  adjustment  routine  (see  Exhibit  2). 

The  computer  program  has  capabilities  to  handle  a  user-supplied 
dependable  capacity  (DCAP)  or  it  can  calculate  dependable  capacity  from  a 
head  and  minimum  discharge  (DC) .  The  program  calculates  the  dependable 
capacity  from  the  discharge  on  the  flow-duration  curve  associated  with  a 
user-supplied  percent  of  time  exceeded  called  SRP,  the  streamflow  relia¬ 
bility  percentage.  The  program  contains  a  default  value  of  85%  for  SRP, 
an  approximation  for  a  15%  forced  outage  percentage  commonly  used  for 
fossil  fueled  generating  plants  which  are  often  the  most  likely  alter¬ 
native  to  hydropower  plants.  The  net  head  is  based  on  the  user-supplied 
discharge  versus  headwater  and  tailwater  relationships. 


9 


The  firm  energy,  FE  expressed  in  megawatt-hours  is  calculated  using 
Equation  9  with  the  maximum  value  of  term  (CAP)j  constrained  to  the 
dependable  capacity. 

(CAP) j  <  DC  (Eq.  11) 

f .  Interruptible  Capacity  and  Energy 


Interruptible  capacity,  IC,  as  used  herein,  is  defined  as  the 
capacity  (KW)  that  the  supplier  can  curtail  at  his  discretion. 

IC  =  CAPDES  -  DC  (Eq.  12) 

Secondary  energy,  SE,  is  simply  the  energy  (MWH)  difference  between  the 
average  energy  and  the  firm  energy. 

SE  =  AE  -  FE  (Eq.  13) 

g .  Potential  Energy  Losses 


The  production  of  average  energy  can  be  limited  by  the 
installed  capacity  of  the  plant,  by  the  discharge  capacity  of  the 
penstock,  or  by  both.  The  program  internally  computes  the  maximum 
potential  energy  by  using  the  power  equation  (Eq.  2)  without  considering 
the  penstock  constraint  (Eq.  3)  or  the  installed  capacity  times  the 
overload  constraint  (Eq.  4).  The  average  energy  based  on  the  power 
equation  (Eq.  2)  with  only  the  installed  capacity  times  the  overload 
constraint  (Eq.  4)  is  calculated  and  subtracted  from  the  maximum 
potential  average  energy  yielding,  AELC,  the  average  energy  lost  due  to 
insufficient  installed  capacity.  Similarly,  the  average  energy 
calculated  with  only  the  penstock  discharge  capacity  as  a  constraint 
(Eq.  3)  is  subtracted  from  the  maximum  potential  average  energy 
yielding,  AELQ,  the  average  energy  lost  due  to  insufficient  penstock 
capacity.  The  constraints  on  the  capacities  used  in  the  calculation  of 
AELC  and  AELQ  are  shown  below  in  Figure  4  for  two  possible  cases. 


Figure  4.  Potential  Energy  Losses 
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1 1 .  Ca 1  cul a t  ion  o_f _  JPower.  Bene fjUt s _ 

HYDUR  can  compute  the  value  cf  power  generated  (power  benefits) 
from  user  supplied  values  (with  several  variations)  or  from  regional  values 
stored  in  the  program.  Capacity  and  energy  benefits  are  specified  as 
functions  of  annual  plant  factor  (the  values  may  be  adjusted  by  ratios). 

The  annual  capacity  benefit  is  computed  as  the  product  of  the  capacity  times 
the  appropriately  interpolated  (based  on  annual  plant  factors)  capacity  value. 
The  annual  energy  benefit,  is  similarly  calculated.  The  total  annual  benefit  is 
the  sum  of  the  two. 

Power  benefits  computed  from  regional  values  stored  in  the  program 
should  be  considered  very  approximate-suitable  for  very  general  preliminary 
studies.  The  benefit  values  are  based  on  regional  studies  performed  by  FKRC 
based  on  1978  data.  Benefit  computations  consider  region  dependable  and 
Lnterrupl ible  capacity  and  firm  and  secondary  energy. 

The  dependable  capacity  benefit  in  dollars  per  year  is  found  by 
multiplying  the  dependable  capacity  by  a  regional  capacity  benefit  factor, 

CB,  which  is  a  function  of  the  plant  factor  (see  Table  la).  The  regions 
applicable  are  defined  in  Figure  5.  The  units  of  CB  are  dollars  per  kilowatt- 
year. 


The  interruptible  capacity  benefit,  in  dollars  per  year  is  computed 
by  multiplying  interruptible  capacity  times  the  product  of  the  capacity 
benefit  and  capacity  benefit  reduction  factor. 

The  energy  benefit  in  dollars  per  year,  is  calculated  by  multiplying 
the  firm  energy  by  the  regional  firm  energy  benefit  value  (see  Table  lb). 

The  secondary  energy  benefit  is  computed  similar  to  the  interruptible  capacity 
benefit . 


The  total  benefit  in  dollars  per  year  is  the  sum  of  the  four  components 
thus  computed. 

i .  Calculation  of  Project  Costs 

HYDUR  can  consider  cost  in  computations  from  user  supplied  values 
or  from  general  relationships  stored  in  the  program.  User  supplied  costs 

are  possible  onLy  for  single  capacity  analyses  -  thus  mav  not  be  used  in  automatic 
sizing  computations.  Program  supplied  cost  functions  are  general  and  are 

therefore  required  for  use  if  optimization  analyses  considering  costs  are  to  hi¬ 
per  formed  . 


Program  supplied  costs  are  taken  directly  from  a  cost  manual  developed 
by  the  Corps  North  Pacific  Division  (NPD(1979)).  The  procedures  were  developed 
for  nation-wide  reconnaissance  level  cost  estimates  of  single-purpose  power 
projects.  The  cost  relationships  were  based  on  empirical  curves  associating 
project  physical  parameters  to  site  component  costs.  All  costs  were  in  July 
1978  dollars.  Costs  generated  from  the  relationships  should  be  considered 
very  preliminary  and  verified  for  reasonableness  subsequent  to  computer  runs. 


11 


Cost  curves  are  included  for  the  powerplant,  embankment,  spillway, 
intake  and  outlet  structures,  waterway,  and  the  reservoir  acquistion  and 
clearing  costs.  Additional  special  cost  items  may  be  included.  Investment 
costs  include  construction  costs  plus  contingency  factor,  engineering  over¬ 
head,  and  interest  during  construction.  A  geographic  adjustment  is  possible. 
Annual  project  costs  are  determined  by  amortizing  these  costs  and  adding 
the  annual  operation,  maintenance,  and  interim  replacement  costs.  A  complete 
description  of  the  cost  estimates  procedures  in  contained  in  the  NPD  (1979) 
document  and  summarized  in  Figure  6. 
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Figure  5.  FERC  Regions  for  Capacity  and  Energy  Benefits 
(from  Hydroelectric  Power  Evaluation,  Federal 
Power  Commission,  Washington,  D.C.,  March  1968 
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COST  ESTIMATE  FORM 

($1,000,000) 

July  1978  Price  Level 
First 

Major  Cost  Items  Cost 

(1)  Powerplant  $ 

(2)  Embankment  (Dams,  Dikes) 

(3)  Spillway 

(4)  Intake  &  Outlet 

(5)  Waterway  (Canal,  Channel,  Conduit) 

(6)  Reservoir 

Investment  Cost  Computations 

(7)  Total  First  Cost  _ 

Subtotal  A  =Z((1)  through  (6))  $ 

(8)  (7)  x  Geographic  Factor  (_.  _)  $ 

(9)  Land  Costs* 

(10)  Subtotal  B  =  ((8)  +  (9))  x  Contingency  (l._  ) 

(11)  Special  Cost  Items** 

(12)  Total  Construction  Cost 
Subtotal  C  =  (10)  +  (11) 


(13) 

(14) 

(15) 

(16) 


Engineering  and  Overhead  Costs 

((12)  x  K  (  .  %) 

engr  —  - 

Total  Project  Cost 
Subtotal  D  =  (12)  +  (13) 


Interest  During  Construction 
((14)  x  IDC  (0. _ )) 

Total  Investment  Cost 
((14)  +  (15)) 


(17) 


Annual  Cost  Computations 
Amortized  Cost 

(Amortization  Factor  (0.  )  x  (16)) 


(18)  Operation  and  Maintenance  Costs 

(19)  Replacement  Costs 

Geographic  Factor  (_.  x  (0.0125)  x  (1) 
x  Contingency  (1.  ) 

(20)  Total  Annual  Cost  =  ((17)  +  (18)  +  (19)) 


Total 

Cost 


$ 


$ 


$ 


$ 


$ 


$ 


*  Geographic  Factor  not  applied  as  Land  requires  no  geographic  adjustment 
for  acquisition. 

**  Geographic  and  Land  Factors  and  Contingencies  not  applied  as  special 

costs  are  user  supplied,  include  concern  for  local  cost  differences,  and 
reflect  level  of  study. 
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Figure  6 


Table  la 


Regional  Energy  Benefit  Valuta  (FERCf  1978) 


FERC  Energy  Benefit  ta  function  of  APF  ($/MWH) 

Region  - - - 


Code 

Region  APP: 

0 

0.10 

0.20 

0 .30 

0.60 

0.50 

0.60 

0.70 

0.80 

0.90 

1.00 

1 

VACAR 

45. J 

45.3 

45.3 

35.7 

35.7 

11.0 

11.0 

4.8 

4.8 

4.8 

4.8 

2 

Southern  Companies  (S.  COMP) 

45.0 

45.0 

45.0 

35.7 

35.7 

9.1 

9.1 

4.8 

4.8 

4.8 

4.3 

3 

ECAR 

38.2 

38.2 

38.8 

23.5 

23.1 

12.7 

12.6 

12.5 

12.4 

12.4 

12.3 

4 

MAIN 

43,9 

43.9 

41.6 

25.4 

23.5 

12.9 

12.4 

12.0 

11.3 

11.6 

11.4 

5 

MARCA 

40.3 

40.3 

37.2 

24.1 

22.6 

10.1 

10.0 

9.9 

9.8 

9.7 

9.7 

6 

WSCC-FW 

33.5 

33.5 

27.4 

24.1 

23.6 

5.8 

6.7 

7.3 

7.0 

8.2 

8.5 

7 

SWPP 

35.2 

35.2 

34.9 

23.3 

22.1 

12.0 

11.9 

11.9 

11.9 

11.8 

11.8 

8 

ERCOT 

29.8 

29.8 

23.8 

22.6 

21.1 

9.4 

9.6 

9.7 

9.8 

9.9 

9.9 

9 

New  England  (NES) 

35.5 

35.5 

30.5 

28.9 

27.1 

1.0 

4.0 

6.0 

7.6 

8.8 

9.8 

10 

New  York  (NYS) 

39.2 

39.2 

39.1 

29.2 

26.9 

10.5 

11.8 

12.8 

13.5 

14.1 

14.5 

11 

PJM 

38.6 

38.6 

36.2 

29.8 

28.1 

11.3 

11.5 

11.6 

11.7 

11.8 

11.9 

12 

CAPCO 

37.6 

37.6 

33.8 

29.8 

26.7 

2,8 

4.5 

5.6 

6.5 

7.2 

7.8 

13 

AEP 

33.2 

33.2 

22.9 

29.8 

25.1 

9.4 

9.5 

9.6 

9.7 

9.7 

9.7 

K 

APS 

45.0 

45.0 

49.1 

31.1 

23.7 

9.4 

9.8 

10.1 

10.3 

10.4 

10.5 

IS 

Northern  California  (N.CA.) 

34.4 

34.4 

35.: 

21.2 

21.6 

11.8 

13.0 

13.8 

34.4 

14.9 

15.3 

16 

Southern  California  (S.CA.) 

33.8 

33.8 

33.6 

21.0 

21.4 

10.0 

11.0 

11.6 

,2.1 

12.5 

12.8 

17 

Pacific  Northwest  (P.N.W.) 

33.9 

31.9 

2b. 7 

21.1 

21.2 

14.1 

13.5 

13.1 

12.7 

12.5 

12.3 

18 

Arizona  (A2) 

34.6 

34.6 

34.6 

21.8 

22  .4 

35.2 

14.7 

14.4 

14.2 

14.0 

13.9 

19 

Southern  Idaho  (S.ID.) 

39.0 

39.0 

43.3 

22.4 

22.4 

9.7 

10.0 

10.2 

10.4 

10.5 

10.6 

20 

Western  Montana  (W.M.) 

39.4 

39.4 

45.2 

22.0 

22.1 

3.4 

3.4 

3.4 

3.4 

3.4 

3.4 

21 

Northern  Nevada  (N.N.) 

38.5 

38.5 

44.8 

21.4 

21.6 

10.7 

11.8 

12.6 

13.1 

13.6 

13.9 

22 

Southern  Nevada  (S.N.) 

36.1 

36.1 

39.5 

21.3 

21.9 

9.3 

9.1 

9.0 

8.8 

6.7 

8.7 

23 

Utah 

34.3 

34.3 

34.0 

20. C 

18.5 

7.7 

7.8 

7.8 

7.8 

7.9 

7.9 

24 

Island  of  Oahu,  Hawaii 

35.5 

35.5 

32.9 

24.2 

25.7 

20.1 

21.2 

22.0 

22.7 

23.1 

23.5 

25 

Island  of  Hawaii,  Hawaii 

7.9 

7.9 

21.4 

25.9 

28.2 

25.7 

26.7 

27.5 

28.0 

28.5 

28.8 

26 

Island  of  Kauai,  Hawaii 

7.6 

7.6 

21.3 

25.9 

28.1 

25.7 

26.7 

27.5 

28.1 

28.5 

28.8 

27 

Island  of  Maul,  Hawaii 

4.4 

4.4 

20.1 

25.4 

28.0 

25.6 

26.8 

27.6 

28.2 

28.7 

29.1 

28 

Island  of  Molakai,  Hawaii 

22.8 

22.8 

32.0 

35.1 

36.6 

37.5 

38.1 

38.6 

38.9 

39.1 

39.3 

29 

Anchorage,  Alaska 

20.6 

20.6 

20.3 

8.7 

11.4 

10.1 

10.9 

11.5 

12.0 

12.4 

12.6 

30 

Fairbanks,  Alaska 

27.0 

27.0 

29.1 

29.8 

30.2 

8.7 

8.9 

9.0 

9.1 

9.2 

9,3 

31 

Valdez,  Alaska 

38.6 

38.6 

38.6 

38.6 

38.6 

38.6 

38.6 

38.6 

38.6 

38.6 

38.6 

32 

Ketchikan,  Alaska 

32.6 

32.6 

32.6 

32.6 

32.6 

32.6 

32.6 

32.6 

32.6 

32.6 

32.6 

Table  lb.  Regional  Capacity  Benefit  Values  (PERC,  1978) 


FERC  Capacity  Benefit  as  function  of  APF  ($/KW) 


Code 

Region  APF: 

0 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

1.00 

1 

VACAR 

25.1 

25.5 

25.9 

40.9 

31.0 

120.1 

121.8 

179.0 

184.3 

189.6 

194.6 

2 

Southern  Companies(S,  COMP)  30.1 

21.2 

12.3 

45.0 

43.1 

109.5 

111.8 

170.7 

175.5 

180.4 

185.3 

3 

ECAR 

31.9 

32.8 

32.8 

66.1 

66.1 

135.2 

135.2 

135.2 

135.2 

135.2 

135.2 

4 

MAIN 

37.2 

33.2 

33.2 

67.1 

67.1 

134.9 

134.9 

134.9 

134.9 

134.9 

134.9 

5 

MARCA 

36.9 

31.5 

31.5 

63.5 

63.5 

135.5 

135.5 

135.5 

135.5 

135.5 

135.5 

6 

WSCC-FW 

40.8 

30.1 

30.1 

68.1 

68.1 

130.8 

130.8 

130.8 

130.8 

130.8 

130.8 

7 

SWPP 

30.8 

30.4 

30.4 

68.9 

68.9 

125.1 

125.1 

125.1 

125.1 

125.1 

125.1 

3 

ERCOT 

39.8 

29.3 

29.3 

65.9 

65.9 

119.0 

119.0 

119.0 

119.0 

119.0 

119.0 

9 

New  England  (NES) 

39.3 

30.5 

?0.5 

70.0 

70.0 

188.1 

188.1 

188.1 

188.1 

188.1 

188.1 

10 

New  York  (NYS) 

33.3 

33.0 

33.0 

75.5 

75.5 

183.7 

183.7 

183.7 

183.7 

183.7 

183.7 

11 

PJM 

32.5 

28.2 

28.2 

64.9 

64.9 

136.0 

136.0 

136.0 

136.0 

136.0 

136.0 

12 

CAPCO 

36.0 

29.3 

29.3 

67.5 

67.5 

180.1 

180.1 

180.1 

180.1 

180.1 

180.1 

13 

AEP 

45.4 

27.3 

27.3 

62.5 

62.5 

110.0 

110.0 

110.0 

110.0 

110.0 

110.0 

14 

APS 

20.4 

27.5 

27.5 

62.5 

62.5 

139.4 

139.4 

139.4 

139.4 

139.4 

139.4 

15 

Northern  Califomia(N.CA.)  35^ 

37.6 

37.6 

69.7 

69.7 

156.5 

156.5 

156.5 

156.5 

156.5 

156.5 

16 

Southe„*n  California(S .CA. )50 .0 

49.6 

49.6 

80.4 

80.4 

164.8 

164.8 

164.8 

164.8 

164.8 

164.8 

17 

Pacific  Notthvest  (P.N.W 

)30.3 

24.7 

24.7 

53.6 

53.6 

121.0 

121.0 

121.0 

121.0 

121.0 

121.0 

18 

Arizona  (AX) 

44.2 

44.3 

44.3 

86.8 

86.8 

224.0 

224.0 

224.0 

224.0 

224.0 

224.0 

19 

Southern  Idaho  (S.ID.) 

21.6 

35.4 

35.4 

71.0 

71.0 

160.1 

160.1 

160.0 

160.1 

160.0 

160.1 

20 

Western  Montana  (W.M.) 

27.1 

37.2 

37.2 

74.3 

74.3 

161.5 

161.5 

161.5 

161.5 

161.5 

161.5 

21 

Northern  Nevada  (N.N.) 

24.1 

35.0 

35.0 

70.1 

70.1 

197.4 

197.4 

197.4 

197.4 

197.4 

197.4 

22 

Southern  Nevada  (S.N.) 

30.5 

36.5 

36.5 

72.6 

72.6 

164.4 

164.4 

164.4 

164.4 

164.4 

164.4 

23 

Utah 

36.8 

^6.2 

36.2 

72.3 

72.3 

162.2 

162.2 

162.2 

162.2 

162.2 

162.2 

24 

Island  of  Oahu,  Hawaii 

49.5 

45.0 

45.0 

75.4 

75.4 

120.6 

120.6 

120.6 

120.6 

120.6 

120.6 

25 

Island  of  Hawaii,  Hawaii 

78.7 

102.3 

102.3 

102.3 

102.3 

169.4 

169.4 

169.4 

169.4 

,  9.4 

169.4 

26 

Island  of  Kauai,  Hawaii 

78.4 

102.3 

102.3 

102.3 

102.3 

169.4 

169.4 

169.4 

169.4 

169.4 

169.4 

27 

Island  of  Maul,  Hawaii 

55.1 

82.7 

82.7 

82.7 

82.7 

171.1 

171.1 

171.1 

171.1 

171.1 

171.1 

28 

Island  of  Molakai,  Hi 

103.6 

119.7 

119.7 

119.7 

119.7 

119.7 

119.7 

119.7 

119.7 

113.7 

119.7 

29 

Anchorage,  Alaska 

30.5 

30.4 

30.4 

48.4 

48.4 

124.5 

124.5 

124.5 

124.5 

124.5 

124.5 

30 

Fairbanks,  Alaska 

35.7 

37.1 

37.1 

37.1 

37.1 

149.5 

149.5 

149.5 

149.5 

149.5 

i.49.5 

31 

Valdez,  Alaska 

109.0 

109.0 

109.0 

109.0 

109.0 

109.0 

109.0 

109.0 

109.0 

109.0 

109.0 

32 

Ketchikan,  Alaska 

49.5 

49.5 

49.5 

49.5 

49.5 

49.5 

49.5 

49.5 

49.5 

49.5 

49.5 
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The  total  annual  cost  developed  from  the  NPD  cost  procedures  is  used 
to  calculate  the  installed  capacity  cost  in  dollars  per  kilowatt-year. 

Likewise  the  average  energy  cost  in  dollars  per  megawatt-hour  is 
computed  by  dividing  the  total  annual  cost  by  the  annual  energy. 

Finally  net  benefits  are  computed  as  total  annual  benefits  minus 
total  annual  cost  and  the  benefit-cost  ratio  determined  as  the  total 
annual  benefit  divided  by  the  total  annual  cost. 


5.  METHODS  OF  OPERATION 

The  HYDUR  program  may  be  operated  in  two  basic  modes.  In  the  first 
mode  the  user  supplies  the  installed  capacity  of  the  plant,  and  the  program 
computes  the  average  annual  energy,  and  if  so  specified,  costs  and  benefits 
as  well.  In  the  second  mode,  the  computer  program  selects  the  installed 
capacity  that  optimizes  a  user  specified  objective  function. 

a •  User  Supplied  Capacity 

In  order  to  invoke  this  option  the  user  simply  provides  the 
capacity  value  CAPDES. 

b.  Optimization  of  Capacity 

When  CAPDES  is  not  specified,  ten  flows  that  correspond  to 
default  intervals  of  1,  5.  10,  20,  40,  60,  90,  95,  99  percent  of  time 
exceeded  on  the  flow-duration  curve  are  used  to  calculate  an  array  of 
capacity  values  nd  their  corresponding  power  parameters.  An  optimi¬ 
zation  scheme  determines  the  installed  capacity  based  upon  a  user 
defined  objective  function. 

When  installed  capacity  currently  exists  at  the  project  all 
calculated  power  parameters  are  incremental  values  (i.e.,  total-existing). 

c .  Data  Entry  Using  the  Alternative  File 

The  user  designates  whether  the  streamflow  or  streamflow 
duration  data  sets  are  to  be  read  from  data  cards  or  from  an 
alternative  file.  The  alternative  file  is  an  existing  data  file  that 
is  created  before  the  HYDUR  program  is  executed.  (See  Exhibit  6.) 

d.  Use  of  the  HYDUR  Program  as  a  Subroutine 

The  HYDUR  program  can  be  used  as  a  separate  subroutine  to 
another  program.  The  requirements  for  doing  so  are  outlined  in  Exhibit  5. 

e .  Overview  of  Program  Options  and  Capab i lity 

I)  Calculates  an  annual  or  seasonal  streamflow  duration 
curve  (SD  card)  from  streamflow  data. 
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2)  Calculates  a  capacity-duration  curve  and  average  energy 
for  any  streamflow  duration  curve. 

-the  installed  capacity  is  either  user  supplied  or 
internally  calculated  by  optimizing  a  user  supplied 
criterion. 

-the  net  head  is  user  supplied  or  both  the  headwater  and 
tailwater  can  be  a  function  of  the  discharge. 

-the  user  supplies  the  minimum  Q  to  turn  the  turbine, 
the  maximum  Q  that  submerges  the  powerhouse,  streamflow 
diversions,  and  any  leakage  or  seepage  through  or  around 
the  dam. 


6.  INPUT 

Input  data  preparation  is  described  in  detail  in  Exhibit  7.  Flow 
data  may  also  be  entered  using  an  alternative  file  as  described  in  Exhibit 
6.  Example  problems  illustrating  input  preparation  are  shown  in  Exhibit  3. 

The  Standard  Input  Format  used  in  the  program  is  ten  8-column  fields 
per  card.  Column  1  and  2  are  used  for  card  identification.  An  integer 
or  alpha  numeric  entry  must  be  right-justified  in  its  field.  A  non-integer 
entry  must  be  right-justified  if  a  decimal  point  is  not  provided.  The 
sequence  of  card  incut  is  generally  insignificant  (except  for  reading  of 
flow  records). 

A  free  format  option  is  available  that  allows  the  user  to  separate 
values  on  an  input  card  with  blanks  or  commas  rather  than  using  the 
standard  10  fields  of  8  columns  each.  This  option  will  facilitate  input 
submitted  through  time  sharing  terminals. 

The  following  cards  are  required  to  perform  a  HYDUR  ananlysis: 

1.  At  least  one  title  card  (T1  through  T4). 

2.  FL  or  FD:  Streamflow  or  Streamflow  Duration  data  cards. 

These  cards  indicate  the  type  of  data  to  be  input  and  its  format. 

3.  PD:  Power  design  data  is  provided  on  this  card.  (Installed 
capacity,  design  head,  efficiency,  etc.) 

4.  EJ :  End  of  Job  Card. 


7  OUTPUT 

Exhibit  4  contains  a  description  of  the  output  data  and  Exhibit  3 
illustrates  output  for  the  example  problems. 
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In  general,  output  from  the  program  is  controlled  through  user 
input  of  the  PS  card.  Detailed  display  of  internal  computer  compu¬ 
tations  may  be  activated  by  proper  coding  of  the  tenth  field  of  this 
card.  (See  Exhibit  7.) 

8.  DEFINITION  OF  TERMS 

Exhibit  7  describes  the  variables  read  as  input  to  the  program, 
while  Exhibit  1  defines  the  important  variables  calculated  within  the 
program. 
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Table  2 


Calculated  Power  Information 


Item 

No. 

Mathematical 

Expression 

Variable 

Name 

Description 

1 

Item 

1 

CAPDES 

Installed  Capacity  (KW) 

2 

Item 

2 

AAE 

Average  Annual  Energy  (MWH) 

3 

Item 

2/(8.76*Item  1) 

APF 

Average  Plant  Factor 

4 

Item 

4 

DC 

Dependable  Capacity  (KW) 

5 

Item 

5 

AFE 

Annual  Firm  Energy  (MWH) 

6 

Item 

1  -  Item  4 

IC 

Interruptible  Capacity  (KW) 

7 

Item 

2  -  Item  5 

ASE 

Annual  Secondary  Energy  (MWH) 

8 

Item 

8 

AEFE 

Annual  Equivalent  Firm  Energy  (MWH) 

-14 

(reserved  for 
future  use) 

15 

Item 

15 

AELQ 

Average  Annual  Energy  Lost  Due  to 
Insufficient  Penstock  Capacity 
(MWH) 

16 

Item 

16 

AELC 

Average  Annual  Energy  Lost  Due  to 
Insufficient  Generating  Capacity 
(MWH) 

17 

Item 

17 

CB 

Dependable  Capacity  Benefit 
($/KW-YR) 

18 

CBR  ■ 

k  Item  17 

ICB 

Interruptible  Capacity  Benefit 
($/KW~YR) 

19 

Item 

19 

EB 

Average  Annual  Energy  Benefit 
($/MWH) 

20 

(Item  4  *  Item  17) 

+  (Item  6  *  Item  18) 

ACB 

Annual  Capacity  Benefit  ($/YR) 

21 

(Item  5  *  Item  19) 

+  Item  7  *  F  (19)) 

AEB 

Annual  Energy  Benefit  ($>/YR) 

18 


Table  2  (Continued) 


Calculated  Power 


Item 

No. 

Mathematical 

Expression 

Variable 

Name 

22 

Item  20  +  Item  21 

TAB 

23 

Item  2 5 /Item  1 

ICC 

24 

Item  25/Item  2 

AAEC 

25 

Item  25 

TAC 

26 

Item  22  -  Item  25 

TANB 

27 

Item  2 2 /Item  25 

BCR 

28-30 

(reserved  for 
future  use) 

Information 

Description 

Total  Annual  Benefit  ($/YR) 
Installed  Capacity  Cost  ($/KW-YR) 
Average  Annual  Energy  Cost  ($/MWH) 
Total  Annual  Cost  ($/YR) 

Total  Annual  Net  Benefit  ($/YR) 
Benefit-to-Cost  Ratio 
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Exhibit  1 


DEFINITIONS  OF  SELECTED  VARIABLES  USED  IN  THE 
HYDUR  PROGRAM 

1.  Variables  contained  in  common  block  POWIN. 


AAEST 

user  supplied  average  annual  energy  estimate  in  MWH. 

This  value  will  override  any  machine-determined 
estimate.  (See  PB  card  in  Exhibit  7.) 

AMOR 

(100)1/ 

Amortization  period  in  years  used  by  the  program  in 
determining  total  annual  costs.  (See  CF  card  in 

Exhibit  7.) 

AVGQ 

machined-determined  estimate  of  average  annual  flow  in 
cfs  or  cms.  This  estimate  is  based  on  integrating  the 
flow-duration  curve  while  accounting  for  any  losses. 

CALCS 

(0) 

indication  code  dealing  with  print  suppression  of  HYDUR 
output.  A  code  of  zero  suppresses  no  output.  A  code  of 
seven  suppresses  output  of  all  calculated  results.  (See 
PS  card  description  in  Exhibit  7.) 

CAPDES 

(0) 

user-supplied  installed  capacity  in  KW.  If  not  defined 
the  program  will  select  a  capacity  value  based  on  some 
objective  function.  (See  PD  card  in  Exhibit  7.) 

CAPYMN  - 
(0) 

user-supplied  minimum  installed  capacity  (KW)  to 
interrogate  when  operating  the  program  in  the 
optimization  mode.  (See  SD  card  in  Exhibit  7.) 

CAPYMX 

(0) 

user-supplied  maximum  installed  capacity  (KW)  to 
interrogate  when  operating  program  in  the  optimization 
mode.  (See  SD  card  in  Exhibit  7.) 

CB 

array  used  to  store  user-supplied  capacity  benefits. 

(See  CB  card  in  Exhibit  7.) 

CBFLAG  - 

logical  variable  indicating  whether  the  program  default 
(true)  or  the  user-supplied  (false)  capacity  benefits 
are  used. 

CBR 

(0.5) 

capacity  benefit  reduction  factor  expressed  as  a  decimal 
fraction.  This  value  is  multiplied  by  the  capacity 
benefit  factor  to  determine  the  value  of  interruptible 
capacity.  (See  PB  card  in  Exhibit  7.) 

1/  Number 

in  parenthesis  is  the  default  value  assigned  to  the 

variable . 
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CEFS 

(0) 


CEMB 

(0) 

CIO  - 

(0) 

CLA 

(0) 

CMIS 

(0) 


CONT 

(0.25) 

COSTR 

(1.0) 


CIVH 

(0) 

CRATIO 

(1.0) 


CRC 

(0) 

CSPW 

(0) 

CWWY 

(0) 

DCAP 

(0) 


DIST 

(0) 


print  suppression  code  dealing  with  cost  estimate  output. 
A  zero  implies  no  print  suppression  while  a  three 
suppresses  all  cost  estimate  output.  (See  PS  card  in 
Exhibit  7.) 

user-supplied  estimate  of  embankment  costs  in  $1,000. 

(See  C2  card  in  Exhibit  7.) 

user-supplied  estimate  of  inlet  and  outlet  costs  in 
$1,000.  (See  C2  card  in  Exhibit  7.) 

user-supplied  estimate  of  land  acquisition  cost  in 
$1,000.  (See  C2  card  in  Exhibit  7.) 

user-supplied  estimate  of  additional  miscellaneous  costs 
in  $1,000  that  are  not  considered  in  the  program's  cost 
routines.  (See  C2  card  of  Exhibit  7.) 

user-supplied  contingency  factor  -to  be  used  in  determin¬ 
ing  total  annual  costs.  (See  CF  card  in  Exhibit  7.) 

adjustment  factor,  expressed  as  a  decimal,  which  is 
applied  to  the  total  annual  cost  calculated  by  the  cost 
routines.  (See  C2  card  in  Exhibit  7.) 

user-supplied  estimate  of  powerhouse  costs  in  $1,000. 

(See  C2  card  in  Exhibit  7.) 

user-defined  adjustment  factor,  expressed  as  a  decimal 
fraction,  which  is  multiplied  to  the  resulting  computed 
total  annual  capacity  benefit.  (See  PB  card  in 
Exhibit  7.) 

user-supplied  estimated  of  reservoir  clearing  costs  in 
$1,000.  (See  C2  card  in  Exhibit  7.) 

user-supplied  estimate  of  spillway  costs  in  $1,000. 

(See  C2  card  in  Exhibit  7.) 

user-supplied  waterway  cost  in  $1,000. 

(See  C2  card  in  Exhibit  7.) 

user-supplied  dependable  capacity  in  KW,  which  will 
override  any  program  estimate.  (See  PB  card  in 
Exhibit  7.) 

user-supplied  dam  length  in  feet  or  meters  required 
to  estimate  embankment  costs  using  the  program's  cost 
routines.  (See  Cl  card  in  Exhibit  7.) 
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DIV 

(0) 


EB 

(0) 

EBFLAG 


ECAP 

(0) 

ECHOS 

(0) 


EFF 

(0.86) 


E RATIO 

(1.0) 


FLOWLO 

(0) 


HEAD 

(0) 

HEIGHT 

(0) 


HMAX 

(100000) 


HMIN 

(0) 


I  COMB 

(0) 


INTBIAS  - 


diversion  in  cfs  or  cms.  Water  diverted  above  the  power 
installation  that  is  not  available  to  produce  power. 

(See  PQ  card  in  Exhibit  7.) 

array  used  to  store  energy  benefits  in  dollars  per  MWH. 
(See  EB  card  in  Exhibit  7.) 

logical  variable  indicating  whether  the  program  default 
(true)  or  the  user-supplied  (false)  energy  benefits  are 
used . 

total  amount  of  existing  installed  capacity  (KW)  at  the 
project.  (See  Cl  card  in  Exhibit  7.) 

print  suppression  code  dealing  with  input  summary 
printed  by  the  program.  A  zero  supresses  no  output 
while  a  one  suppresses  all  input  summary.  (See  PS  card 
in  Exhibit  7.) 

combined  efficiency  of  the  turbine  and  generator  units 
expressed  as  a  decimal  fraction.  EFF  will  be  overridden 
if  TE  cards  are  input.  (See  PD  card  in  Exhibit  7.) 

user-supplied  adjustment  factor,  expressed  as  a  decimal 
fraction,  which  is  multiplied  to  the  resulting  computed 
total  annual  energy  benefit.  (See  PB  card  in  Exhibit  7.) 

indicator  which  specifies  whether  or  not  a  valid 
ctreamflow  value  has  been  input.  (See  FL  card  in 
Exhibit  7.) 

gross  operating  power  head  in  feet  or  meters. 

(See  PD  card  in  Exhibit  7.) 

user-supplied  dam  height  in  feet  or  meters  required  to 
estimate  embankment  costs  using  the  programs's  cost 
routine.  (See  Cl  card  in  Exhibit  7.) 

user-supplied  power  head  in  feet  or  meters  which  is  the 
upper  operating  limit  for  the  installed  unit(s).  (See 
PQ  card  in  Exhibit  7.) 

user-supplied  power  head  in  feet  or  meters  which  is  the 
lowest  operating  limit  for  the  installed  units(s). 

(See  PQ  card  in  Exhibit  7.) 

combination  code  used  in  the  defining  of  an  objective 
function  for  selecting  an  optimal  capacity  size.  (See 
0C  card  in  Exhibit  7.) 

initial  pointer  to  the  streamflow  values  used  in 
developing  seasonal  flow-duration  curves.  (See  SD  card 
in  Exhibit  7.) 
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IPOWANL  - 


IPRNTS 

(0) 

IPROJ 

(0) 


IREG 

(0) 

ISEASN  - 

(1) 

IVAR1 

(26) 


IVAR2 

(0) 


JSTATE 

(0) 


KS 

(1) 

LENREC  - 
(80) 

LSTRD 

NCLUDE 

(ALL) 

NJOB 

NINYR 


NOP 


flag  which  indicates  whether  a  complete  power  analysis 
or  only  a  flow-duration  curve  be  constructed.  (See  FL 
card  in  Exhibit  7.) 

print  suppression  code  used  in  seasonal  power  analysis 
runs  only.  (See  PS  card  in  Exhibit  7.) 

code  used  to  delete  consideration  of  certain  primary 
costs  components  calculated  by  the  cost  routine.  (See 
Cl  card  in  Exhibit  7.) 

region  code  to  select  the  proper  regional  benefit  curves 
in  the  program.  (See  PB  card  in  Exhibit  7.) 

contains  the  current  season  number  (1-12)  during  seasonal 
power  analyses  runs  only. 

item  description  code  which  corresponds  to  the  first  data 
item  used  in  defining  the  objective  function.  (See  0C 
card  in  Exhibit  7.) 

item  description  code  which  corresponds  to  the  second 
data  item  used  in  defining  the  objective  function.  (See 
OC  card  in  Exhibit  7.) 

state  code  used  by  the  cost  routine  to  determine  the 
proper  geographical  adjustment  to  apply  to  the 
calculated  primary  cost  components.  (See  C2  card  in 
Exhibit  7.) 

valley  slope  code  used  in  determining  embankment  costs. 
(See  Cl  card  in  Exhibit  7.) 

record  length  of  the  alternative  file.  (See  FL  Card 
in  Exhibit  7. ) 

temporarily  not  used. 

number  of  consecutive  flow  values  used  in  developing  a 
flow-duration  curve.  (See  SD  card  in  Exhibit  7.) 

current  job  number  of  the  entire  computer  run. 

number  of  flow  values  comprising  a  year.  (See  SD  card 
in  Exhibit  7 . ) 

number  of  points  comprising  the  OP  array. 
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NPQ 

(70) 


an  integer  variable  describing  the  number  of 
coordinates  used  to  define  the  streamflow  duration  curve. 


NSEASKP  -  number  of  flow  values  to  skip  between  seasonal  power 

(0)  runs.  (See  SD  card  in  Exhibit  7.) 

NSEASN  -  number  of  seasons  to  perform  in  power  analyses. 

(1)  (See  SD  card  in  Exhibit  7.) 

NT  -  number  of  points  comprising  the  tailwater-discharge 

relationship.  (See  TV  card  in  Exhibit  7.) 

OP  -  array  of  percent  exceedance  values  used  in  the  capacity 

selection  mode  of  the  program.  (See  OP  card  in 
Exhibit  7.) 

OPER  -  operational  mode  indicator  which  is  used  in  the 

(L)  determination  of  annual  operation  costs.  (See  Cl  card 

in  Exhibit  7 . ) 

OPERND  -  indicator  which  determines  if  the  objective  function 

(MAX)  used  in  selecting  capacity  should  be  maximized  or 

minimized.  (See  0C  card  in  Exhibit  7.) 

OVLOAD  -  installed  capacity  overload  factor.  (See  PD  card  in 

(1.15)  Exhibit  7.) 

PE  -  array  used  to  store  percent  exceedence  of  the 

streamf low-duration  curve. 

PEAKF  -  peaking  factor  used  to  analyze  the  marketability  of  a 

(1)  project.  (See  PD  card  in  Exhibit  7.) 

PLOTS  -  print  suppression  code  dealing  with  plotted  output  of  the 
(0)  program.  (See  PS  card  in  Exhibit  7.) 

PSR  -  power  storage  ratio  used  in  analyzing  storage  project. 

(0)  (See  PD  card  in  Exhibit  7.) 

PTC  -  time  of  project  construction  in  years.  Used  to  estimate 

(2)  interest  during  construction.  (See  CF  card  in 

Exhibit  7.) 

QDES  -  total  plant  design  capacity  in  cfs  or  cms.  (See  PD  card 

(0)  in  Exhibit  7. ) 

0FACT  -  decimal  fraction  that  is  used  to  adjust  the  streamflow 

(1.0)  ordinates  of  the  streamf low-duration  curve.  (See  PQ 

card  in  Exhibit  7.) 

QMIN  -  minimum  operating  discharge  limit  for  the  installed 

(0)  unit(s).  (See  PQ  card  in  Exhibit  7.) 
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QQ 


QSUB 

(0) 

RAAE 

(1.0) 


RATE 

(.06875) 

RE  PL 
(.0125) 

RESA 

(0) 


SEAEND 

(FALSE) 

SEBR 

(-50) 


SRP 

(-85) 

T 


TE 


TERM 

(-999) 

TH 


TITLE 


TQ 


array  containing  the  streamflow  value  of  the  streamflow- 
duration  curve.  (See  QQ  ca^d  in  Exhibit  7.) 

discharge  value  where  tailwater  effects  preclude  power 
operation.  (See  PQ  card  in  Exhibit  7.) 

ratio,  expressed  as  a  decimal  fraction,  which  is  applied 
to  the  machine-determined  average  annual  energy  value. 
(See  PB  card  in  Exhibit  7.) 

discount  rate  used  in  amortizing  total  investment  costs. 
(See  CF  card  in  Exhibit  7.) 

replace  cost  factor,  expressed  as  a  decimal  fractional 
part  of  powerhouse  costs.  (See  CF  card  in  Exhibit  7.) 

user-supplied  reservoir  area  in  acres  or  square 
kilometers  used  :'n  the  deteimination  of  reservoir 
clearing  costs.  (See  Cl  t  ari  in. Exhibit  7.) 

logical  variable  indicating  whether  a  seasonal  power 
analysis  is  complete.  (.TRUE,  implies  yes). 

energy  benefit  reduction  factor,  expressed  as  a 
decimal  fraction,  which  is  applied  to  the  firm  energy 
benefit  factor  to  determine  secondary  energy  benefits. 
(See  PB  card  in  Exhibit  7.) 

streamflow  reliability  factor  used  in  the  determination 
of  dependable  capacity.  (See  PB  card  in  Exhibit  7.) 

logical  array  that  Indicates  whether  any  of  the  four 
title  cards  (Tl,  T2,  T3  or  T4)  in  the  TITLE  (20,4)  array 
have  been  previously  defined.  When  T(i)  is  true  the 
program  will  print  the  Tj  card  contained  in  the  TITLE 
(20, i)  array.  T(i)  should  be  specified  as  false  when 
the  Tj  card  is  not  defined. 

array  containing  efficiencies  that  wil]  be  related  to 
power  discharge.  (See  TE  card  in  Exhibit  7.) 

termination  code  used  in  reading  streamflow  values 
from  QQ  cards.  (See  FM  card  in  Exhibit  7.) 

array  containing  headwater  elevations  that  will  be 
related  to  reservoir  inflows.  (See  TH  card  in 
Exhibit  7.) 

array  containing  title  information  supplied  by  the  user 
on  Tl  through  T4  cards. 

array  containing  flow  values  that  will  be  related  to 
information  contained  in  the  TE,  TH,  and  TO  arrays. 
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TRACE  -  cumulative  trace  code  used  for  displaying  internal 
(0)  computer  computations. 

TURB  -  turbine-type  indicator.  (See  Cl  card  in  Exhibit  7.) 

(0) 

TW  -  array  containing  tailwater  elevations  that  will  be 

related  to  turbine  discharges.  (See  TW  card  in 
Exhibit  7.) 

UAPF  -  plant  factor  related  to  firm  energy  delivery.  (See 

(1.0)  PD  card  in  Exhibit  7.) 

WYL  -  user-supplied  waterway  length  used  to  determine  waterway 

(0)  costs.  (See  Cl  card  in  Exhibit  7.) 

WYQ  -  user-supplied  waterway  design  flow  used  to  determine 

(0)  waterway  costs.  (See  Cl  card  in  Exhibit  7.) 

XTRP1  -  extrapolation  variable  used  in  powerhouse  cost  routine. 

XTRP2  -  extrapolation  variable  used  in  powerhouse  cost  routine. 

XUNITS  -  number  of  units  comprising  total  installed  capacity  of 

(0)  the  plant.  (See  PD  card  in  Exhibit  7.) 

2.  Variables  used  in  common  block  POWOUT.  These  variables  are 
returned  to  this  common  block  for  each  capacity  size  tested. 

AAE  -  average  annual  energy  in  MWH. 

AAEC  -  average  annual  energy  cost  in  dollars  per  MWH. 

ACB  -  annual  capacity  benefit  in  dollars  per  year. 

AEB  -  annual  energy  benefit  in  dollars  per  year. 

AEFF  -  efficiency  of  plant. 

AELC  -  average  annual  energy  that  is  lost  due  to  insufficient 

generating  capability  in  MWH. 

AELQ  -  average  annual  energy  that  is  lost  due  to  insufficient 

penstock  capacity  in  MWH. 

AFE  -  annual  firm  energy  in  MWH. 

AHEAD  -  average  power  head  in  feet  or  meters. 

AHEADW  -  average  headwater  elevation  in  feet  or  meters. 

ADP  -  annual  plant  factor  expressed  as  a  decimal  fraction. 
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AQ 


average  annual  inflow  in  cfs  or  cms. 


AQGEN  -  average  annual  flow  used  for  generation  in  cfs  or  cms. 

ASE  -  average  annual  secondary  energy  in  MWH. 

ATAILW  -  average  tailwater  elevation  in  feet  or  meters. 

BCR  -  benefit  to  cost  ratio. 

CAPCTY  -  installed  capacity  determined  by  the  program  in  KW. 

CBX  -  dependable  capacity  benefit  in  i/KW-YR. 

DC  -  dependable  capacity  in  KW. 

DHEAD  -  design  head  in  feet  or  meters  associated  with  the 
installed  capacity. 

EBX  -  firm  energy  benefit  in  $/MWH. 

IC  -  interruptible  capacity  in  KW. 

ICB  -  interruptible  capacity  benefit  in  $/KW-YR. 

ICC  -  installed  capacity  cost  in  $/KW-YR. 

QDES  -  design  discharge  in  cfs  or  cms  associated  with  the 

Installed  capacity. 

TAB  -  total  annual  benefits  in  $/YR. 

T/C  -  total  annual  costs  in  $/YR. 

TANB  -  total  annual  net  benefits  in  $/YR. 

EQF  -  equivalent  firm  energy  in  MWH. 
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Discharge 


EXHIBIT  2 


ADJUSTMENT  OF  FLOW-DURATION  CURVE 
FOR  STORAGE  EFFECTS 

A)  DEVELOPMENT  OF  METHODOLOGY 

The  analysis  of  hydropower  storage  projects  has  traditionally  been 
performed  by  use  of  sequential  reseivoir  routing  techniques.  The  use 
of  flow-duration  techniques  has  been  traditionally  applied  only  to 
run-of-river  type  projects. 

While  individual  power  storage  projects  should  be  analyze.-  by 
detailed  sequential  routings  when  sufficient  funds  and  detail  are 
available,  the  flow-duration  technique  (as  modified  herein)  can  be 
made  to  somewhat  approximate  the  results  of  a  sequential  routing  by 
modifying  the  flow-duration  curve  to  represent  outflow  conditions. 

A  storage  project,  in  general,  accumulates  excessive  inflows  for 
future  use  during  low-flow  periods,  thereby  transforming  the 
inflow-duration  curve,  based  on  inflows  into  the  project,  into  a 
flatter  outflow-duration  curve,  reflecting  the  operation  and  effect  of 
the  project's  storage  as  depicted  below: 


FLOW-DURATION  CURVE 
Based  on  INFLOWS 


Area  =  A0 


+ 


STORAGE 

EFFECT 


H - — - 

0  %  Exceedence 


Figure  ) 
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Superimposing  the  previous  curves  results  in: 


Figure  2 


where  thp  area  (A0)  under  the  original  flow-duration  curve  is 
preserved  and  the  modified  flow-duration  curve  passes  through  points 
Cj  and  Cp,  where 

Cjl  represents  a  discharge  corresponding  to  100  percent  exceedence; 

C2  represents  a  point  of  intersection  between  the  two  curves. 

The  following  analytical  technique  was  developed  to  determine  the 
transformed  shape  of  the  flow-duration  curve  that  satisfies  the  above 
conditions.  The  criteria  used  in  making  a  selection  of  and  C2 
will  be  discussed  later. 

Restated,  the  problem  requires  the  development  of  an  algorithm 
that  will  generate  a  modified  flow-duration  curve  and  meet  the 
following  conditions: 

1)  the  value  of  the  function  (flow-duration  ordinate)  at  100 
percent  exceedence  must  be  Cj ; 

2)  the  value  of  the  function  (flow-duration  ordinate)  at  some 
percent  exceedance  p  must  be  C2,  where  0  <p  <1.0; 

3)  the  area  under  the  modified  flow-duration  curve  must  equal 
the  area  under  the  original  flow-duration  curve  (Aq). 
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The  first  assumption  made  in  the  technique  is  to  decide  the 
mathematical  form  of  the  function  to  be  used.  Since  flow-duration 
curves  resemble  decaying  power  functions  the  following  generalized 
function  was  assumed: 


Q(x) 


AD 


(B(l-x) ) 


Eq.  1 


where  A  =  multiplying  coefficient; 

D  =  a  base  value; 

B  =  a  power  coefficient; 

and  where  x  =  the  percent  exceedence  corresponding  to  discharge  (Q)> 
where  x  is  expressed  as  a  fraction  between  0  and  1. 


After  initial  testing,  it  was  found  that  the  assumed  power 
function  could  only  satisfy  any  two  of  the  three  conditions  required. 
Therefore,  the  analytical  form  was  modified  to: 


Q(x)  =  A  10(B(1-x)} . Eq.  : 

where  A  =  multiplying  coefficient; 

B  =  a  power  coefficient; 

and  where  x  =  the  percent  exceedence  correspondence  to  discharge  (Q), 
where  x  is  expressed  as  a  fraction  between  0  and  1. 

Eq.  2  can  now  be  used  to  define  a  relationship  which  passes 
through  points  and  C2  designated  as  Ql(x).  The  same 

mathematical  form  can  be  used  to  develop  another  relationship 
designated  as  Q2(x),  which  passes  through  C2  in  such  a  manner  as  to 
ensure  the  area  under  both  curves  is  Aq. 

The  necessary  conditions  can  be  restructured  as  follows: 

1)  Define  a  function  (Ql(x))  such  that  Ql(l.O) 

=  Cj  and  Ql(p)  =  C2', 

2)  Determine  the  area  under  Ql(x)  from  x  =  p  to  x  “  1,0  and 
call  the  area  A^; 

3)  Define  another  function  (Q2(x))  such  that  Q2(p)  =  C2 
and  the  area  under  the  curve  from  x  =  0  to  x  =  p 
(designated  as  A2)  is  equal  to  the  original  area  A0 
minus  Aj . 
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The  following  illustrates  the  above  requirements: 


Figure  3 


The  coefficients  A  and  B  for  function  (Ql(x))  can  be  determined  by 
satisfying  condition  (1): 

A  =  Ci . Eq.  3 

B  =  100(02/0!)/  (1  -  p) . Eq.  4 

Integrating  Ql(x)  from  x  =  p  to  x  =  1.0  yields  the  following  for 
condition  (2): 


Ax=  (C1/(2.3B))(10(B(1"p))-1)  . Eq.  5 

To  this  point  in  the  technique,  a  relationship  for  the  modified 
flow-duration  curve  from  x  =  p  to  x  =  1.0  has  been  developed  and  the 
area  under  this  segment  of  the  curve  has  been  determined.  The 
function  (Q2(x))  must  satisfy  the  remaining  condition  (3): 
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Integrating  Q2(x)  from  x  =  0  to  x  =  p,  and  noting  that  the 
function  (Q2(x))  must  pass  through  Cj  at  x  =  p  yields: 

A2  =  A0  -  Ai . Eq.  6 

2.3A2/C2  =  [10(BP)-1]  /B . Eq.  7a 

There  is  no  direct  solution  for  B  in  Equation  7a;  however,  a  value 
for  B  can  be  obtained  by  use  of  Newton's  first  order  approximation. 

This  trial  and  error  solution  is  conditional  upon  C2  being  within  an 
acceptable  range  (i.e.,  there  exists  some  values  MN  and  MX,  such  that 
MN <  C2 <  MX  is  a  condition  necessary  for  solution  of  B). 

A  requirement  frr  the  use  of  the  first-order  approximation 
technique  is  an  estimate  of  the  first  derivative  of  Equation  7a. 

Closer  inspection  of  Equation  7a,  reveals  the  expression  can  be 
simplified  since  A2  and  C2  have  been  previously  determined: 

F(Bp)  =  ((10<BP)  -  1)/B)  -  K . Eq.  7b 


where  K  =  2.3A2/C2 

Substituting  the  dummy  variable  u  =  Bp  yields: 

F(u)  =  p((10u  -  l)/u)  -  K . Eq.  8 

Taking  the  first  derivative  of  F(u)  results  in  the  following: 

dF (u)  =  p  (2.3ulOu  -  10u  +  1)/  u2 . Eq.  9 

du 

Of  interest  is  where  the  first  derivative  equals  zero  (corresponds  to 
the  minimum  value  of  u) .  Solving  Equation  9  for  zero  results  in  u 
approaching  zero  in  the  limit,  since  at  u  =  0,  Equation  9  is  undefined. 

An  algorithm  can  now  be  devised  to  generate  the  modified  flow- 
duration  curve  as  outlined  below,  given  the  area  A0  and  parameters 
Ci  and  C2-  (Knowing  C2  implies  knowing  p). 

1)  Solve  the  coefficients  A  and  B  in  power  function  (Ql(x))  by 
using  estimates  Ci  and  C2  in  Equations  3  and  4; 

2)  Determine  the  area  under  function  (Ql(x))  using  Equation  5 
and  call  result  A2; 

3)  Determine  the  remaining  required  area  needed  under  function 
(02(x))  from  Equation  6  and  call  result  A2; 
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4)  Calculate  the  K  component  of  Equation  7b  using  results  from 
step  3  and  the  estimate  of  C2 ; 

5)  Test  to  see  if  the  minimum  value  for  u,  say  u  =  0.01,  results 
in  a  value  less  than  K.  If  not,  make  new  estimate  of  C2 

and  return  to  Step  1; 

6)  Use  Equation  9  to  make  new  estimate  of  u  and  continue  process 
until  the  estimate  of  u  results  in  Equation  8  equaling  zero. 

B)  PARAMETER  DETERMINATION 

Attention  is  now  focused  on  making  estimates  for  parameters  C^, 

C2  and  p.  The  value  of  Aq  is  a  constant  and  represents  the  area 
underneath  the  original  flow-duration  curve,  which  is  easily 
determined  by  integration. 

The  value  of  C2  has  been  shown  to  be  critical  in  a  final 
feasible  solution  of  the  problem.  This  parameter  will  be 
systematically  set  equal  to  Aq,  which  generally  corresponds  to 
percent  exeedences  ranging  in  value  from  15  to  35  percent.  The 
selection  for  C2  will  automatically  determine  the  value  of  p  since 
C2  and  p  are  functionally  related  through  use  of  the  flow-duration 
relationship.  Therefore,  a  means  of  selecting  a  value  for  Cj  is  the 
only  obstacle  remaining  in  defining  a  plausible  solution. 

The  value  of  Cj  is  dependent  on  the  storage  capability  of  the 
site  being  analyzed.  Accordingly,  it  seems  reasonable  to  assume  that 
Cj  can  be  estimated  by  considering  the  base  flow  component  of  the 
flow  regime  and  the  minimum  flow  contribution  due  to  reservoir 
regulation  during  adverse  flow  conditions  as  follows: 

Cx  *  QMC  +  QMSC . Eq.  10 


where : 

Cj  =  the  minimum  flow  value  on  the  synthetic  flow-duration 
curve  corresponding  to  100  percent  exceedence; 

QMC  =  the  minimum  flow  value  on  the  original  flow-duration 
curve  without  regard  to  storage  effects  (100  percent 
exceedence  value); 

QMSC  =  the  minimum  flow  contribution  attributed  to  reservoir 
operation  under  critical  low  inflow  conditions. 
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Critical  low  flow  conditions  occur  whenever,  over  a  sustained 
period  of  time,  a  reservoir  is  regulated  to  release  additional  flow  in 
excess  of  upstream  inflows  as  a  means  of  satisfying  designed  project 
purposes.  With  regard  to  hydropower,  this  operational  policy,  if 
continued,  can  actually  exaggerate  the  situation  since  depletion  of 
power  storage  reduces  the  effective  headwater  and  correspondingly  the 
operating  power  head;  requiring  a  continually  increasing  amount  of 
flow  to  sustain  energy  requirements.  The  period  of  maximum  drawdown 
can  be  defined  as  the  period  of  time  which  begins  with  full  power 
storage  and  ends  when  the  power  storage  remaining  is  at  a  minimum.  By 
definition,  the  period  of  maximum  drawdown  will  then  contain  the  most 
adverse  streamflow  conditions  and  will  require  the  maximum  withdrawal 
of  water  from  the  power  storage.  An  estimate  of  QMSC  can  now  be 
approximated  by  determining  the  depletion  rate  occurring  throughout 
this  period . 

As  an  initial  step,  the  power  storage  can  be  converted  from  units 
of  volume,  typically  in  acre-feet,  to  units  of  flow  rate  (cfs).  To 
perform  this  conversion,  a  time  period,  say  one  year,  must  be 
selected.  The  resulting  value  expresses  the  power  storage  potential 
as  the  average  amount  of  flow  that  can  be  extracted  from  an  initially 
full  power  pool  throughout  a  period  of  one  year. 

However,  as  defined  above,  the  period  of  maximum  drawdown  is 
unconstrained  with  regard  to  the  length  of  time  required  to  complete 
the  process,  and  actual  reservoir  operations  have  demonstrated  this 
period  of  time  varies  from  a  few  weeks  to  several  years  in  length. 
Accordingly,  the  initial  depletion  rate  (assuming  a  one-year  length  in 
the  period  of  maximum  drawdown)  must  be  adjusted  by  a  factor  to 
reflect  the  project's  actual  length  in  time  to  minimum  pool  level  as 
shown  below: 

QMSC  =  PS  *  ACF  *  ADJF . Eq.  11 

where 

PS  =  power  storage  expressed  as  a  volume  (acre-feet); 

ACF  =  a  conversion  factor  (0.00138)  which  when  multipled  by  (PS) 
expresses  the  amount  of  power  storages  in  terms  of  an 
average  annual  flow  rate  (CFS-YR); 

ADJF  =  adjustment  factor  applied  to  the  power  storage  to  correct 
for  variation  in  the  length  of  the  period  of  maximum 
drawdown. 

From  Equation  11,  one  can  conclude  that  a  length  of  the  period  of 
maximum  drawdown  exceeding  one  year  requires  the  adjustment  factor 
(ADJF)  to  be  less  than  one,  since  a  factor  failing  this  requirement 
will  cause  QMSC  to  be  overestimated.  Conversely,  a  length  in  the 
period  of  maximum  drawdown  less  than  unity  requires  ADJF  to  exceed 
unity.  In  effect,  ADJF  can  be  alternatively  defined  as  the  reciprocal 
of  the  length  in  the  period  of  maximum  drawdown,  when  time  is  measured 
in  years. 
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Several  attempts  at  establishing  a  relationship  to  determine  a 
value  for  ADJF  were  performed.  The  regression  equation  finally 
selected  is  based  upon  313  existing  and  proposed  hydro  sites 
throughout  the  United  States  and  is  shown  below: 

ADJF  =  0.65  +  1.113*  L0G(1/PSR) . Eq.  12 

where  PSR  represents  a  project's  power  storage  to  mean  annual  flow  ratio. 

Statistically,  the  above  equation  resulted  in  a  R-squared  value  of 
0.49  and  a  standard  error  of  0.45.  A  plot  of  this  relationship  can  be 
seen  in  Figure  4. 

The  PSR  is  a  dimensionless  parameter  which  expresses  the  relative 
size  of  power  storage  to  average  annual  inflow  and  is  determined  by 
converting  power  storage  to  an  average  one  year  flow  rate,  as 
previously  suggested,  and  then  dividing  the  result  by  the  project's 
expected  average  annual  inflow. 

A  relatively  large  PSR,  say  greater  than  1.0,  indicates  sufficient 
storage  capacity  so  that,  on  the  average,  there  exists  the  capability 
of  extending  the  period  of  maximum  drawdown  during  sustained  periods 
of  low  inflow.  As  the  PSR  falls  below  1.0,  this  capability  to 
attenuate  diversity  between  inflow  availability  and  project  demands 
decreases,  causing  the  average  length  of  the  period  of  maximum 
drawdown  to  decrease,  accordingly.  Another  observation  in  Figure  4, 
substantiating  this  conclusion,  is  that  as  the  PSR  falls  below  1.0, 
considerable  increase  in  scatter  of  the  data  occurs,  implying  that  the 
decreasing  storage  capability  is  becoming  relatively  less  effective  in 
controlling  the  diversity  between  inflow  supply  and  energy  demand.  A 
direct  determination  of  parameter  is  now  possible  by  successive 
use  of  Equations  12,  11,  and  10;  allowing  for  a  plausible  solution  to 
systematically  producing  a  synthetic  flow-duration  relationship  for 
projects  exhibiting  power  storage. 

The  substitution  of  the  relationship  for  the  original  flow- 
duration  curve  will  substantially  improve  any  estimate  of  average 
energy  and  will  additionally  enable  an  approximation  of  dependable 
capacity  to  be  performed  when  used  in  a  nonsequential  power  potential 
analysis. 

Dependable  capacity  can  be  defined  as  the  capacity  which  under  the 
most  adverse  flow  conditions  on  record,  can  be  relied  upon  to  carry 
system  load,  provide  dependable  reserve  capacity,  and  meet  firm  power 
obligations,  taking  into  account  seasonal  variations  and  other 
characteristics  of  the  load  to  be  supplied.  The  association  to  the 
"most  adverse  flow  conditions  on  record,"  in  the  definition  of 
dependable  capacity  strongly  supports  the  notion  that  parameter  Cj 
might  be  valuable  as  an  indicator  in  approximating  this  capacity 
value.  This  assumption  was  tested  and  it  was  found  that  dependable 
capacity  could  be  estimated  by  using  the  power  equation  as  follows: 
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2  8  II  ADJF  =  0.85  +  1.1 1  LOG  (1/PSR) 
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FIGURE  4  PLOT  OF  REGRESSION  ANALYSIS  RESULT 


DCAP  =  C(Ci/PF)He 


Eq.  13 


where : 

DCAP  =  dependable  capacity  in  kilowatts; 

C  =  .084603  conversion  factor  which  expresses  power  in 
kilowatts; 

Cj  =  the  minimum  flow  parameter  as  previously  defined; 

PF  =  the  average  annual  plant  factor  relating  dependable 

capacity  to  its  firm  energy  requirement; 

H  =  net  power  head  in  feet; 

e  =  overall  efficiency  (assumed  equal  to  a  constant  0.86) 

The  quantity  (Cj/PF)  represents  the  expected  minimum  flow  adjusted 
for  average  time  of  hydropower  plant  operation. 

Equation  13  was  employed  with  data  from  all  113  projects  used  to 
develop  the  parameter  ADJF  and  the  resulting  estimates  of  dependable 
capacity  were  compared  to  corresponding  dependable  capacity  values 
estimated  from  sequential  routing  techniques. 

This  comparison,  depicted  in  Figure  5,  resulted  in  an  R-squared 
value  of  0.985  and  a  standard  error  of  10,700  kilowatts.  Comparison 
of  empirical  to  sequentially  determined  capacities  varied  over  a  range 
of  100  kilowatts  to  650  megawatts.  A  departure  in  plotting  position 
above  an  imaginary  45°  line  represents  an  underestimate  of  dependable 
capacity  determined  empirical  as  compared  to  dependable  capacity 
estimated  using  sequential  routing  techniques.  Conversely,  a 
departure  below  this  line  represents  an  exaggeration  of  dependable 
capacity.  Maximum  departure  about  this  line,  in  terms  of  percent 
difference,  occurs  for  small  installations  (i.e.,  project  having 
installed  capacities  less  than  2  megawatts).  Small  installations  are 
generally  associated  with  projects  possessing  limited  storage  capacity 
and  correspondingly  small  power  storage  to  mean  annual  flow  ratios 
(PSR).  Since  the  regression  equation  used  to  estimate  ADJF  (Equation 
12)  was  incorporated  in  the  empirical  determination  of  dependable 
capacity,  the  problem  of  increased  scatter  associated  with  small  PSR's 
is  most  probably  the  underlying  influence  causing  these  maximum 
departures  to  occur  in  this  capacity  range. 
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Example  Problem  Illustrating  The  Use  of 
Synthetic  Flow-Duration  Curve  Techniques 


A  :ume  the  following  hydrologic  data  can  be  developed  for  a  site: 

1)  power  storage  (PS)  =  42,783  acre-feet 

2)  net  power  head  (H)  =  148  feet 

3)  dependable  capacity  plant  factor  APF  =  0.6017 

4)  estimate  of  dependable  capacity  from  sequential  analysis 
=  1;049  Kw 

5)  original  flow-duration  curve  (see  Figure  6) 

A)  Estimate  of  Dependable  Capacity 

Step  1  -  Integrating  the  flow-duration  curve  yields  an 
average  annual  flow  (AAQ)  of  95.22  cfs 

Step  2  -  Determine  power  storage  to  mean  flow  ratio  (PSR) 

HSR  =  PS  *  0.00138/AAQ 

PSR  =  42783  *  .00138/95.22  =  0.62 

Step  3  -  Use  PSR  and  regression  equation  to  find 
ajustment  factor  ADJF 

ADJF  *  0.65  +  1.113  L0G(1/PSR)  =  0.88 

Step  4  -  Calculate  minimum  contribution  of  flow  due  to 
storage  regulation  (QMSC) 

QMSC  =  PS*. 00138  *  ADJF 

QMSC  =  42783  *  .00138  *  .88  =  52.02  cfs 

Step  5  -  Find  total  minimum  flow  (Cj) 

Cj.  =  QMC  +  QMSC 

C2  =  1.05  +  52.02  =  53.07  cfs 
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Example  Problem  Illustrating  The  Use  of 
Synthetic  Flow-Duration  Curve  Techniques 
(continued) 


Step  6  -  Determine  dependable  capacity  DCAP 

DCAP  =  (Cj/ADF)  *  H  *  C  *  e 

DCAP  =  (53.07/0.6017)  *  148  *  0.084603  *  0.86 
=  951  KW 

Note:  Ratio  of  dependable  from  sequential  and  algorithm  =  1049/951  =  1.10 

B)  Synthetic,  flow-duration  curve  algorithm 

Step  1  -  Select  cross-over  point  between  curves  (C2) 
at  percent  exceedence  corresponding  to  AAQ 

C2  «  AAQ  =  95.22  cfs 

Note:  C2  corresponds  to  percent  exceedence  (p).  Approximately 
30  percent  (see  Figure  6) 

Step  2  -  Determine  decaying  power  function  Ql(x)  from 
Ql(.30)  to  Ql(1.0) 

A  =  Cl  =  53.07 

B  =  100(02/0! )/(l-p) 

B  =  LOG  (95.22/53.07)(l-0.3) 

B  =  0.36 

therefore  Ql(x)  =  53.07*10**(0.36(l-x) ) 

Step  3  -  Determine  area  under  Ql(x)  from  x  =  0.30 
to  x  =  1.0 

A1  =  (Cj./ (2.3  *  B))  *  (10**(B(l-p) )-l) 

Ax  =  (53.07/(2.3  *  0.36)  *  (10**(0. 36( .7) )-l) 

A!  =  50.59  cfs 
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Example  Problem  Illustrating  The  Use  of 
Synthetic  Flow-Duration  Curve  Techniques 
(continued) 

Step  4  -  Determine  area  under  Q2(x)  from  x  =  0  to  x 
x  =  .30 

A2  =  AAQ  -  Ax  -  95.22  -  50.59  =  44.63  cfs 

Step  5  -  Determine  decaying  power  function  Q2(X)  from 
x  =  0  to  x  =  0.30 

K  =  2.3  *  A2/C2 

K  =  2.3(44.63)/95.22  =  1.08 

using  Newton's  1st  order  approximation  solve  for 
F(u)  =  0 


where 


F(u)  =  p((10u-l)/u)  -  K 
-  F'(u)  =  p(2.3ulOU  -  10U  +  l)/u2 
and  p  =  0.30 


1st  Guess:  try  u  =  1.0 

F(1.0)  =  .3(10-1)/1  -  1.08  =  1.62 
since  F(1.0)  i  0,  try  again 
F'  (1.0)  =  0.3(2.3(10)  -  10  +  1)/1  =  4.20 
2nd  Guess  u  =  u  -F(u)/F'  (u)  =  1  -  1.62/4.20  ■  0.61 
F(0. 61)  =  0.3(10*61  -  l)/0. 61  -  1.08  =  0.432 
since  F(0.61)  ^  0,  try  again 

F'  (0.61)=0.3(2.3(.61)  lO’61  -  10'61  +1)/(.61)2  =  2.130 

F'  (u)  implies  the  first  derivative  of  function  F(u)  with  respect 
to  u. 
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Example  Problem  Illustrating  The  Use  of 
Synthetic  Flow-Duration  Curve  Tecuniques 
(continued) 

3rd  Guess  u  =  0.61  -  0.432/2.130  =  0.41 

F(0.41)  =  0.3  (10*41  -l)/0. 41  -  1.08  =  0.069 
since  F(0.41)  40  try  again 

F'  (0.41)  =  0.3  (2.3( .41)  10*41  -10*41  +  1)/(.41)2  =  1.523 

4th  Guess  u  =  0.41  -  .069/1.523  =  0.36 

F(0. 36)  =  0.3(10’36  -l)/.36  -  1.08  =  -0.004 

since  F(0.36)  is  %  0,  iterative  process  completed.  Recalling 
u  =  Bp  and  p  =  0.30,  then; 

B  =  0.36/. 3  =  1.20 

By  definition,  at  x  =  p,  Q2(x)  =  C2.  Therefore,  by 
substitution  of  these  parameters  in  function  (Q2(x))  yields: 

A  =  C2/10(B  "Bp)  =  95.22/10*84  =  13.76 

Parameter  A  and  B  for  Q2(x)  are  now  determined,  therefore, 

Q2(x)  =  13.76  *  io(1*2°(1-x>) 

Use  of  Ql(x)  and  Q2(x)  will  duplicate  synthetic  curve  depicted  in 
Figure  6. 
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GRAPHICAL  REPRESENTATION  OF  EXAMPLE 


EXHIBIT  3 


TEST  PROBLEMS 


The  teat  problems  of  this  exhibit  are  included  to  illustrate  selected 
analytical  capabilities,  input  requirements,  and  output  format  of  the  HYDUR 
program.  The  problems  are  also  intended  for  use  in  verification  of  distributed 
program  code.  Four  problems  are  included  to  demonstrate  a  range  of  potential 
applications  of  the  HYDUR  program.  Specific  problem  statements,  description 
of  input  requirements,  and  results  are  provided  for  each  problem. 

1.  TEST  PROBLEM  -  STTE  ASSESSMENT  FOR  ENERGY  CAPACITY 

This  test  problem  presents  a  basic  application  of  the  HYDUR  program. 

The  problem  depicts  a  preliminary  assessment  of  the  feasibility  of  installing 
100,000  kilowatt  (PI). 4)  generating  capacity  at  an  existing  reservoir  designed 
and  operated  for  irrigation  purposes.  Two  units  (PD. 8)  of  50,000  kilowatts 
each  are  analyzed.  The  flow  duration  curve  is  input  directly  in  standard 
format  of  10F8.0  as  specified  on  the  FD  card.  The  last  value  of  the  flow- 
duration  data  is  followed  by  a  -999  which  indicates  the  end  of  the  data  set. 

The  following  assumptions  were  made  with  respect  to  this  test  problem. 

1)  The  maximum  penstock  discharge  capacity  is  input  as  7710  cfs  (PD.l). 

2)  The  powerplant  efficiency  is  estimated  to  be  .86  (PD. 2). 

3)  An  overload  factor  of  1.2  (PD. 5)  was  input. 

4)  An  upstream  diversion  of  47.35  cfs  (PQ.3)  exists  that  includes  a 
water  supply  requirement  and  net  evaporation.  This  loss  will  be 
subtracted  from  available  power  flow  values  but  will  not  affect 
power  operating  head. 

5)  A  diversion  loss  (leakage)  of  10  cfs  was  input  (PQ.4)  to  account 
for  flow  losses  that  would  not  contribute  to  power  production 
(through  the  penstocks)  but  WOuld  affect  tailwater  elevations. 

NOTE:  For  this  problem  with  a  constant  input  head  of  175  feet  tail- 
water  elevations  will  not  fluctuate. 

6)  An  adjustment  factor  of  .958  (PQ.5)  was  input  to  account  for  the 
difference  in  discharge  value  (drainage  area  ratios)  between  the 
downstream  gage  and  the  damsite.  This  value  is  multiplied  times 
each  streamflow  value. 

7)  A  set  of  TQ,  TH,  and  TW  cards  were  used  to  define  tailwater  conditions. 

The  results  of  the  analysis  indicate  that  132,977  MwH  of  average  annual 
energy  efuld  be  generated  by  the  two  50,000  KW  units.  The  facility  would  use 
a  calculated  average  outflow  of  1200  cfs  for  power  generation. 
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TABLE  OF  DISCHARGE  (CFS/CMS)  VS.  EFFICIENCY,  HEADWATER  AND  TAILWATER  (FT/KT) 

TQ  CARD  0.00  3000.00  9000.00  105000.00  289000.00  578000.00  1156000.00  2312000.00  4624000.00  9248000.00 
TE  CARD  0.860  0.860  0.860  0.860  0.860  0.860  0.860  0.860  0.860  0.860 
TH  CARD  577.50  577.50  577.50  577.50  577.50  577.50  577.50  577.50  577.50  577.50 
TW  CARD  400.00  404.00  408.00  440.00  457.00  457.00  457.00  457.00  457.00  457.00 
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PRELIMINARY  ANALYSIS  OF  POTENTIAL  HYDROPOWER  SITE 

TWO  UNITS  WITH  TOTAL  CAPACITY  OF  100000KW  CAPACITY 
STREAMFLOWS  IN  STANDARD  FORMAT,  USGS  GAGE  7339000,  1930-1967 
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PERCENT  OF  TIME  EXCEEDED 


PRELIMINARY  ANALYSIS  OF  POTENTIAL  HYDROPOWER  SITE 

TWO  UNITS  WITH  TOTAL  CAPACITY  OF  100000KW  CAPACITY 
,Y  STREAMFLOWS  IN  STANDARD  FORMAT,  USGS  GAGE  7339000,  1930-1967 
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2.  TEST  2  -  FEASIBILITY  ASSESSMENT  OF  INSTALLING  HYDROPOWER  AT  EXISTING 

STORAGE  PROJECT. 

This  test  problem  assesses  the  potential  maximum  hydropower  capacity  at 
an  existing  storage  project  which  presently  has  no  generating  facilities. 

The  optimization  capability  of  the  program  is  used  in  this  analysis. 

The  ananlysis  uses  an  input  flow-duration  curve  (FD  Card).  The  maximum 
penstock  discharge  capacity  PD.l  was  left  blank  since  the  maximum  installed 
capacity  (0C.1)  was  to  be  determined,  CAPDES  (PD. 4  =  blank).  The  design 
plant  efficiency  (PD. 2)  was  input  as  .86,  overload  factor  (PD. 5)  as  1.2,  the 
annual  plant  factor  (UAPF)  as  .137,  and  the  power  storage  adjustment  factor 
(PD. 7)  as  .49.  Two  units  were  assumed  to  be  installed  (PD. 8).  The  flow 
diversion  (PQ.3)  of  47.35  cfs  was  assumed,  with  50  cfs  estimated  to  lost 
(PD. 4)  due  to  leakage  and  a  fish  ladder  operation.  Tailwater  fluctuation 
affects  on  head  were  defined  by  TQ,  TH,  and  TW  cards. 

The  results  indicate  a  maximum  installed  capacity  of  31,947  KW  yielding 
an  average  annual  energy  of  123,345  MWH.  The  natural  and  adjusted  flow  duration 
curves  are  defined  in  the  output,  which  is  included,  along  with  the  input  on 
the  following  pages. 
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CARD  0.00  3000.00  9000.00  10500.00  289000.00  578000.00  1156000.00  2312000.00  4624000.00  9248000.00 
CARD  0.860  0.860  0.860  0.860  0.860  0.860  0.860  0.860  0.860  0.860 
CARD  577.50  577.50  577.50  577.50  577.50  577.50  577.50  577.50  577.50  577.50 
CARD  400.00  404.00  408.00  440.00  457.00  457.00  457.00  457.00  457.00  457.00 
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PRINT  SUPPRESSION 


ANALYSIS  OF  POTENTIAL  HYDROPOWER  STORAGE  PROJ1 
DURATION  CURVE  ADJUSTED  FOR  STORAGE  AFFECT! 
iPACITIES  CALCULATED  FOR  SPECIFIED  FREQUENCIES 
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PERCENT  OF  TIME  EXCEEDED 
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FLOW  DURATION  CURVE  ADJUSTED  FOR  POWER  STORAGE  EFFECTS 

ADJUSTED  FLOW  DURATION  CURVE  —  FLOWS  IN  CFS/CMS  AND  PERCENT  OF  TIME  EXCEEDED  ORDINATES  AS  DECIMAL  FRACTIONS 
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ANALYSIS  OF  POTENTIAL  HYDROPOWER  STORAGE  PROJEi 
DURATION  CURVE  ADJUSTED  FOR  STORAGE  AFFECTS 
CAPACITIES  CALCULATED  FOR  SPECIFIED  FREQUENCIES 
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ANALYSIS  OF  POTENTIAL  HYDROPOWER  STORAGE  PROJECT 
DURATION  CURVE  ADJUSTED  FOR  STORAGE  AFFECTS 
CAPACITIES  CALCULATED  FOR  SPECIFIED  FREQUENCIES 
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3. 


TEST  PROBLEM  3  -  ECONOMIC  FEASIBILITY  OF  ADDING  POWER  TO  MONTECELLO  DAM 


This  test  problem  evaluates  the  economic  feasibility  of  adding  hydro¬ 
electric  power  to  an  existing  project  which  has  no  present  generating  capacity. 

The  objective  is  to  determine  automatically  (using  the  optimization  capability 
of  the  program)  the  total  potential  capacity  which  minimizes  the  cost  of  pro¬ 
ducing  energy  expressed  in  terms  of  dollars  per  megawatt  hours  (MWH) . 

The  major  input  specifications  for  the  problem  include:  the  flow  duration 
relationship,  power  facility  design  and  flow  information,  tailwater  specifications, 
and  optimization-criteria.  The  site  has  a  maximum  penstock  discharge  capacity  of 
1124  cfs  (PD . 1) ,  an  upstream  diversion  of  99.4  cfs  (PQ.3),  and  a  minimum  flow  for 
power  operation  of  75  cfs  (PQ.2).  An  average  head  of  400  feet  is  specified  in 
field  PD. 3,  with  the  tailwater  -  discharge  relationship  provided  on  the  TQ  and 
TW  cards.  Headwater  elevations  are  determined  by  the  program  by  adding  400  feet 
(PD. 3)  to  the  TW  card  values  instead  of  inputting  headwater  elevations  on  the  TH 
card.  Thirteen  years  of  monthly  streamflow  data  developed  at  the  site  (PQ.5  =  1.0) 
and  input  in  a  user-supplied  format  of  10F8.0  (FM.3,  FM.4-10)  are  provided.  The 
FERC  regionalized  benefits  are  overriden  by  user-supplied  benefits  provided  on  the 
EB  and  CB  cards.  All  costs  of  the  existing  project  are  surpressed  except  for  the 
cost  of  the  powerhouse  (Cl. 10). 

The  optimization  criteria  specified  was  to  minimize  the  cost  of  producing 
energy  (OC.l  =  24  and  0C.4  =  MIN).  A  22.2  percent  time  of  exceeded  flow  values 
satisfies  this  objective  with  an  average  annual  energy  cost  of  $21.01  per  MWH. 

The  determined  capacity  to  be  installed  is  13,939  KW  which  is  capable  of  generating 
60,754  MWH  of  energy  annually.  Note  that  the  default  option  of  maximizing  net 
benefits  would  have  altered  the  installed  capacity  selected. 

The  input  and  output  for  Test  3  are  shown  on  the  following  pages. 
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UNRECOGNIZED  VALUE  FOR  OPER  ON  Cl  CARD  —  L  AS! 
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MONTE  CELLO  DAM 

POWER  ANALYSIS  AT  AN  EXISTING  DAM 
MONTHLY  FLOWS  IN  USER  FORMAT,  USGS  GAGE  11454000,  1959-1971 
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PERCENT  OF  TIME  EXCEEDED 
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CITY 


EXTRAPOLATION  OF  COST  DATA  REQUIRED 

INTAKE/OUTLET  COSTS  INHIBITED  DUE  TO  CAPACITY  LT  25 , 000 
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POWER  ANALYSIS  AT  AN  EXISTING 
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4. 


TEST  PROBLEM  4  -  THE  DALLES  RUN-OF-R1VER  PROJECT 


This  test  problem  analyzes  the  economic  feasibility  of  adding  capacity  ti 
an  existing  1,806,800  KW  (Cl. 4)  run~of-river  project  (PD. 7=0).  Both  the  pen¬ 
stock  capacity  and  the  plant  capacity  are  to  be  optimized  so  the  PD.l  and  PD. 4 
fields  are  left  blank.  A  negative  net  evaporation  of  -12  cfs  (precipitation 
gain)  is  entered  in  field  PQ.3.  The  net  power  head  versus  discharge  relation¬ 
ship  is  provided  on  the  TH  and  TQ  cards.  Since  power  head  and  not  headwater 
elevations  were  use  on  the  TH  cards,  tailwater  elevations  of  zero  must  be  entered 
on  the  TW  card.  An  annual  plant  factor  for  firm  energy  delivery  is  assumed  to 
he  1.0  (PD. 6). 

The  streamflow-durat ion  curve,  developed  from  streamflows  at  site  (PQ. 5=1.0) 
is  entered  via  cards  (FD.1=CARD)  in  standard  format  (FD.2=STD).  Default  energy 
and  capacity  benefits  for  the  FERC  region  code  17  (PB.l)  are  used.  Costs  are 
based  on  information  piovided  on  the  Cl  and  C2  cards  using  the  North  Pacific 
Division,  (NPD)  Corps  of  Engineers,  cost  estimating  procedures.  The  calculated 
total  annual  cost  is  adjusted  by  1.3  (C2.10).  The  OP  card  specifies  that  the 
optimization  summarizes  tables  for  the  5,  5.5,  6.  8.  10,  12,  14,  16,  18,  and  20 
percent  of  time  ordinates  (existing  capacity  corresponds  to  approximately  20  per¬ 
cent  time  of  exceeded  o-dinate).  The  PS  card  suppresses  all  plots  (PS.l)  and  the 
tabulation  of  the  streamf low-duration  curve  (PS. 4). 

The  optimization  table  for  various  percent  of  time  exceeded  displays  two 
columns  exhibiting  a  negative  incremental  capacity  (corresponding  to  18  and  20 
percent  exceedence).  The  total  annual  cost,  9999999999,  is  used  by  the  cost 
estimating  routines  whenever  an  estimate  of  powerhouse  costs  is  not  possible. 

The  two  parameters  needed  for  this  cost  component  estimate  are  the  unit  capacity 
and  corresponding  design  head.  (Refer  to  the  NPD  cost  estimating  manual  for 
examples  of  powerhouse  cost  curves). 

The  additional  installed  capacity  determineu  from  the  optimization  analysis 
is  144,061  KW  which  corresponds  to  an  exceedence  value  of  14.3  percent  of  the 
time.  The  total  capacity  of  the  facility  would  be  1,950,861  KW.  The  incremental 
additional  capacity  is  economically  feasible  to  install  assuming  given  benefits 
and  costs. 

The  input  and  output  for  Test  4  are  shown  on  the  following  pages. 


TEST  4 
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UNRECOGNIZED  VALUE  FOR  OPER  ON  Cl 
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EXHIBIT  4 


DESCRIPTION  OF  PROGRAM  OUTPUT 

The  output  from  the  program  falls  into  the  following  categories: 

(1)  banner  page;  (2)  flow-duration  input;  (3)  summary  of  the  input 
variables  and  selected  job  options;  (4)  a  plot  of  the  unadjusted  flow- 
duration  curve;  (5)  ordinates  and  plot  of  the  adjusted  streamflow-duration 
curve;  (6)  cost  estimate  from  the  existing  capacity;  (7)  plot  of  the 
capacity-duration  curve  for  the  existing  capacity;  (8)  extrapolation 
warning  messages;  (9)  summary  of  optimization  calculations;  (10)  cost 
estimate  form  for  the  total  potential  capacity;  (11)  summary  of  the 
power  analysis  for  total  potential,  existing,  and  incremental 
capacity;  (12)  and  a  plot  of  the  capacity-duration  curve  for  the  total 
potential  capacity.  A  detailed  description  of  the  items  that  appear 
in  the  output  is  given  below, 

(1)  Banner  Fage.  The  banner  page  contains  the  program  title.  On 
the  Corps  of  Engineers  library  versions  at  the  Lawrence 
Berkeley  Laboratory  (LBL)  and  Boeing  Computer  System  (3CS)  a 
general  and  HYDUR  message  are  printed.  Respectively,  these 
provide  general  information  about  the  Corps  computer  library 
and  specific  comments  about  the  HYDUR  program. 

(2)  Flow-Duration  Input.  This  page  provides  a  copy  of  the  input 
cards  for  the  entire  run. 

(3)  Summary  of  Input  Variables  and  Selected  Job  Options.  All  the 
input  variables  and  selected  job  options  used  in  the 
indicated  job  number  are  summarized.  When  several  jobs  are 
stacked,  this  summary  will  indicate  what  options  are  in 
effect  for  the  current  job. 

(4)  Plot  of  the  Flow-Duration  Curve.  The  unadjusted  streamflo'w 
duration  plot  is  provided  with  notes  indicating  QLOSS,  QMIN, 
and  DIV.  The  streamflows  are  plotted  on  a  logarithmic  scale 
versus  an  linear  percent  of  time  scale. 

(5)  Adjusted  Flow-Duration  Curve  Ordinates.  The  unadjusted 
streamflow-duration  ordinates  are  multiplied  by  QFACT  and  DIV 
is  subtracted  yielding  the  adjusted  streamflow  ordinates  that 
are  printed.  If  PD. 7  is  positive  the  curve  is  also  adjusted 
for  power  storage  effects  (see  Exhibit  2). 

(6)  Cost  Estimate  Form  for  Existing  Capacity.  The  cost  estimate  form 
for  the  existing  capacity  is  printed  when  the  value  for  ECAP 

is  non-zero  (Cl. 4).  This  form  is  based  on  the  cost 
estimating  procedures  developed  by  the  North  Pacific  Division 
(NPD).  Warning  messages  are  printed  above  the  form  when  the 
cost  functions  are  extrapolated. 
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(7)  Plot  of  Capacity-Duration  Curve  for  Existing  Capacity.  A 
plot  of  the  capacity-duration  curve  is  printed  when  ECAP  is 
non-zero  (Cl. 4).  The  capacity  values,  which  are  limited  to 
ECAP  times  the  overload,  are  plotted  on  a  logarithmic  scale. 

(8)  Extrapolation  Warning  Messages.  Warning  messages  are 
printed  when  extrapolation  occurs  in  the  cost  routines  that 
are  called  during  the  construction  of  the  optimization 
calculations.  The  program  internally  cycles  20  times 
through  the  optimization  table  so  some  messages  may  be 
repeated  up  to  20  times. 

(9)  Summary  of  Optimization  Table.  The  optimization  table 

summarizes  the  thirty  items  described  in  Table  2  of  the  users 
manual.  These  items  are  printed  for  the  capacities 
corresponding  to  the  percent  of  time  exceeded  ordinates 
requested  on  the  OP  card.  When  ECAP,  the  existing  capacity, 
is  specified  (Cl. 4),  all  items  displayed  are  incremental 
values  (total  potential  minus  the  existing  potential 
values) .  The  lower  portion  of  the  table  summarizes  the 
average  operating  characteristics  in  the  power  analysis 
which  include  the  inflow  available  for  generation,  the 
outflow  used  for  power  generation,  the  efficiency,  net  head, 
headwater  elevation,  and  tailwater  elevation.  In  addition 
the  design  penstock  capacity  is  printed.  This  value  is 
either  the  value  supplied  by  the  user,  or  when  not  supplied 
it  is  computed  by  the  program  to  provide  just  enough  flow  to 
allow  generation  of  the  installed  capacity.  The  design  head 
is  calculated  from  the  installed  capacity  and  the  design 
flow.  At  the  bottom  of  the  table  the  optimum  percent  of 
time  exceeded  is  printed.  This  is  based  on  the  optimization 
criterion  provided  on  the  0C  card. 

(10)  Cost  Estimate  Form  for  the  Total  Potential  Capacity.  The 
cost  estimate  for  the  total  potential  capacity  is  based 
either  on  CAPDES  (PD. 4)  when  it  is  supplied,  or  on  the 
optimized  installed  capacity.  Any  cost  extrapolation 
warning  messages  are  printed  above  the  cost  estimate  form. 

(11)  Summary  of  the  Power  Analysis.  The  summary  of  the  power 
analysis  tabulates  che  thirty  items  in  Table  2  of  the  users 
manual  for  the  potential,  existing  and  incremental 
capacities.  The  potential  capacity  column  is  based  on 
either  the  optimized  total  potential  or  CAPDES  when  it  is 
provided  (PD. 4).  When  there  is  no  existing  capacity 

the  potential  and  incremental  capacity  columns  will  be  the 
same.  When  the  potential  capacity  is  greater  than  a 
non-zero  existing  capacity  (Cl. 4),  the  incremental  capacity 
column  is  based  on  the  difference.  If  the  optimized  total 
potential,  which  is  calculated  internally,  is  less  than  the 
existing  capacity,  then  the  printed  total  potential  column 
Ip  set  equal  to  the  existing  capacity  column. 
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(12)  Plot  of  the  Capacity-Duration  Curve  for  Total  Potential 

Capacity.  A  plot  of  the  capacity-duration  curve  is  printed 
for  the  total  potential.  CAPDES  on  the  plot  is  based  on 
either  the  value  supplied  in  field  (PD. 4)  or  is  the  optimized 
installed  capacity.  (Note  that  the  asterisks  on  the  plot 
override  the  alphanumeric  characters). 
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EXHIBIT  5 


USE  OF  HYDUR  AS  A  SUBROUTINE 

The  HYDUR  PROGRAM  can  be  used  as  a  subroutine  to  an  existing 
program.  For  the  purposes  of  this  discussion,  the  users  existing 
program  will  be  referred  to  as  EXIST.  To  accomplish  this  task  the 
main  program  of  the  HYDUR  PROGRAM  called  DRIVER  is  deleted  and  two 
labelled  common  blocks,  type  declarations,  and  a  call  statement  are 
added  to  the  program  EXIST.  The  first  common  block  called  P0WIN  as 
shown  below  supplies  the  hydropower  routines  with  all  the  input 
variables  that  would  normally  be  read  on  cards  if  HYDUR  were  run  as  a 
stand-alone  program. 

COMMON/POWIN/AAEST,  AMOR,  AVGQ,  CALCS,  CAPDES ,  CB ,  CBFLAG ,  CBR,  CEFS, 
.CEMB,  CIO,  CLA,  CMIS,  CONT,  COSTR,  CPWH,  CRATIO,  CRC,  SEBR,  CSPW, 
.SEAEND,  REPL,  CWWY,  DCAP ,  DIST,  DIV,  EB,  EBFLAG,  ECAP,  ECHOS,  EFF, 
.ERATIO,  HEAD,  HEIGHT,  ICOMB,  IPROJ,  IREG,  IVAR1,  IVAR2 ,  J STATE ,  KS, 
.NJOB,  NOP,  NPQ,  NT,  OP,  OPER,  OPERND,  OVLOAD,  PE,  PLOTS,  PSR,  PTC, 
•QDES,  QFACT ,  QLOSS,  QMIN,  QQ,  QSUB,  RAAE,  RATE,  RESA,  SRP,  IPRTS, 
.FLOWLO,  T,  TE,  TH,  TITLE,  TQ,  TRACE,  TURB,  TW,  UAPF,  WYL,  WYQ,  XUNITS, 
.PEAKF ,  HMIN,  HMAX,  NINYR,  INTBIAS ,  NCLUDE,  NSEASN,  NSEASKP ,  ISEASN, 
.IPOWANL,  LENREC,  LSTRD,  CAPYMN,  NRDS,  NPRO,  CAPYMX,  TERMCD ,  XTRPl, 
.XTRP2 

INTEGER  CALCS,  CEFS,  ECHOS,  OPER,  OPERND,  PLOTS,  TITLE(20,4),  TRACE, 
•TURB,  XTRPl,  XTRP2 

REAL  CB(ll) ,  EB(ll),  OP(ll),  PE(70),  QQ(70), 

•TE(10),  TH(10) ,  TQ(10) ,  TW(10) 

LOGICAL  CBFLAG,  EBFLAG,  T(4).  SEAEND 

All  of  these  variables  are  described  in  Exhibit  1. 

The  second  common  block  P0W0UT  contains  all  the  important  variables 
calculated  in  HYDUR.  The  definition  of  these  variables  is  contained 
in  Exhibit  1. 

COMMON/ POWOUT/  AAE,  AAEC ,  ACB ,  AEB ,  AEFF,  AELC,  AELQ,  AFE,  AHEAD, 

. AHEADW ,  APF ,  AQ,  AQGEN,  ASE,  ATAILW,  BCR,  CAPCTY,  CBX,  DC,  DHEAD , 

•EBX,  IC,  ICB ,  ICC,  QDSIGN,  TAB,  TAC,  TANB ,  EQF 

REAL  IC,  ICB,  ICC 

The  call  statement  is  as  follows: 

CALL  HYDUR  (SWITCH) 

When  SWITCH,  a  logical  variable,  is  ’TRUE*  the  variables  contained 
within  P0WIN  must  be  defined  in  EXIST  prior  to  the  call  statement. 

If  SWITCH  is  'FALSE'  then  the  data  cards  normally  read  as  input  to  the 
HYDUR  program  must  be  supplied.  The  location  of  the  HYDUR  data 
relative  to  the  data  read  by  EXIST  will  be  dependent  on  where  the  call 
statement  to  HYDUR  was  inserted  in  program  EXIST. 
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DATA  ENTRY  USING  THE  ALTERNATIVE  FILE 

Sets  of  streamflow  data  or  streamf low-duration  coordinates  can  be 
read  into  the  program  from  an  alternative  file.  The  first  field  on 
the  FL  and  FD  cards  (see  Exhibit  7)  indicates  hether  the  data  sets 
are  to  be  read  by  cards  (using  option  CARD)  or  from  the  alternative 
file  (using  option  FILE). 

The  alternative  file  is  an  external  file  that  must  be  created 
prior  to  the  execution  of  the  HYDUR  PROGRAM.  The  alternative  file 
consists  of  all  the  data  sets  requested  by  the  FILE  opticn  on  the  FL 
and  FD  cards.  The  data  sets  must  be  located  on  the  alternative  file 
in  the  same  order  as  requested  by  the  program.  Each  line  of  the 
alternative  file  consists  of  a  card  image  which  cannot  exceed  150 
characters  in  length. 

The  program  reads  data  from  the  alternative  file  on  Unit  3.  If 
the  user's  alternative  file  is  named  XYZ,  the  user  must  associate  this 
external  (physical)  file  XYZ  with  the  internal  (logical)  file  named 
Unit  3.  This  task  is  accomplished  through  job  control  language  (JCL) 
which  varies  from  computer  to  computer.  Examples  1,  2  and  3  describe 
how  this  is  accomplished  for  three  computer  installations. 

Figure  J  illustrates  an  example  mutliple  job  run  which  requests 
three  data  sets  from  the  alternative  file.  Job  1  on  figure  1  reads 
daily  streamflow  data  (FL.l  specifies  DAILY)  from  the  alternative  file 
(FL.2  specifies  FILE).  After  the  program  reads  the  FL  card  it  commences 
to  read  data  from  the  alternative  file  until  a  different  streamgage 
number  is  encountered  in  the  next  set  of  header  cards.  In  this 
example,  this  is  accomplished  by  inserting  a  blank  H  and  N  card  after 
the  last  daily  streamflow  card.  The  second  job  reads  flow-duration 
data  in  STD  format  from  the  alternative  file  until  the  terminator  code 
of  -999  is  encountered.  Streamflow  data  are  read  from  cards  in  the 
third  job  which  demonstrates  that  jobs  reading  data  sets  from  cards 
can  be  intermingled  with  jobs  reading  data  sets  from  the  alternative 
file.  The  fourth  job  reads  streamflow  data  in  STD  format  until  the 
-999  is  encountered. 

After  each  data  set  is  read  from  the  alternative  file,  the  file 
remains  positioned  to  read  the  next  data  requested  by  the  FL  or  FD 
card  specifying  the  FILE  option.  End-of-record  or  end-of-file  marks 
are  not  inserted  between  data  sets. 
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Figure  l.  Use  of  the  Alternative  File  to  Read  Flow  Data 


Example  1:  JCL  for  CDC  Computers  at  Boeing  Computer  Services 

Methods  A  and  B  yield  equivalent  results.  Methods  assumed  file  (XYZ) 
is  in  users  library. 

Line  Method  A  Method  B 


1 

2 

3 

4 

5 

6 


GET,  TAPE  3  =  XYZ. 

GET,  HYDUR/UN  =  CECELB. 
HYDUR. 

/*EOR  (end-of-record  card) 
data  cards 

/*EOF  (end-of-file  card) 


GET,  XYZ. 

GET,  HYDUR/UN  =  CECELB. 
HYDUR,,,,,, XYZ. 

/*EOR  (end-of-record  card) 
data  cards 

/*£0F  (end-of-file  card) 


Method  A 

Line  1  creates  a  local  copy  of  alternative  file  XYZ  called  TAPE  3. 

This  matches  the  program’s  internal  name  for  this  file. 

Method  B 

Line  1  creates  a  local  copy  of  alternative  file  XYZ.  In  line  3,  file 
substitution  is  used  to  create  the  association  between  internal  file 
TAPE  3  and  external  file  XYZ.  (TAPE  3  is  the  sixth  file  named  on  the  CDC 
"program"  card  in  HYDUR). 
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Example  2:  JCL  for  the  CDC  computers  at  Lawrence  Berkeley  Laboratory 
(LBL)  in  Berkeley,  California. 

Methods  A  and  B  yield  equivalent  results. 


Line 


Method  A 


Method  B 


1 

2 

3 

4 

5 

6 


FETCHPS , yourli b, TAPE3 , XYZ . 
FETCHPS ,HECLIB , HYDUR ,HYDUR. 
HYDUR. 

7/8/9  (end-o f-record  card) 
data  cards 

6/ 7/8/9  (end-of-file  card) 


FETCHPS ,yourl ib, XYZ ,XYZ . 
FETCHPS ,HECLIB , HYDUR , HYDUR. 
HYDUR,,,,,, XYZ. 

7/8/9  (end-of-record  card) 
data  cards 

6/7/ 8/ 9  (end-of-file  card) 


Method  A 

Line  1  creates  a  local  copy  of  alternative  file  XYZ  called  TAPE  3 
which  is  located  on  a  user  defined  library  (yourlib)  in  a  subset, 
called  XYZ.  The  name  TAPE  3  matches  the  program's  internal  name  for 
this  file. 

Method  B 


Line  1  creates  a  local  copy  of  alternative  file  XYZ.  In  line  3,  file 
substitution  is  used  to  create  the  association  between  internal  file 
TAPE  3  and  external  file  XYZ.  (TAPE  3  is  the  sixth  file  named  on  the 
CDC  "program"  card  in  HYDUR.) 


6 / 7 / 8 / 9  -  multipunched 
7/8/9  -  multipunched 


Exhib 
4  of 


*h  m 


Example  3:  JCL  for  IBM  S360/370  computers 

Only  the  basic  JCL  structure  is  shown.  The  user  must  supply  other 
file-dependent  and  installation-dependent  parameters. 

//STEPA  EXEC  PGM=HYDUR 

//FT01F001  DD  DCB= (LRECL=  ) 

//FT02F001  DD  DSN=  , DCB=(RECFM=F,LRECL=80,  ) 

//FT06F001  DD  SYS0UT=A 
//FT05F001  DD  * 
data  cards 

/* 

Brief  Description  of  FTO  files: 

FT01F001  Scratch  file 
If  this  is  a  formated  write: 

DCB=(RECFM=1 ,LRECL=m,  ) 

where  m  is  the  record  size  in  bytes. 

If  this  is  an  unformated  write  and  the  record  size  can  vary 
(within  or  between  runs)  then: 

DCB=( RECFM=V , LRECL=n ,  ) 

where  n=maximum  bytes  per  record  plus  eight. 

FT02F001  Alternative  data  file 

FT06F001  Printed  output 

FT05F001  In-job  card  input 
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HYDUR  INPUT  DESCRIPTION 


Exhibit  7 


Input  Description 


September  1982 


This  exhibit  contains  a  detailed  description  of  each  input  card  used 
by  computer  program  "HYDUR" .  Many  of  the  cards  shown  can  be  omitted  if 
certain  options  are  not  required.  The  summary  of  input  cards  at  the  end 
of  this  exhibit  shows  the  sequential  arrangement  of  cards  and  the  location 
of  variables  on  the  cards.  The  location  of  variables  for  each  input  card 
is  shown  by  field  number. 


The  standard  format  of  input  cards  is  ten  fields  of  eight  columns 
each  except  for  Field  1 .  Variables  in  Field  1  occupy  card  columns  3-8 
because  card  columns  1  and  2  (referred  to  as  Field  0)  are  reserved  for  the 
required  identification  characters.  Exceptions  are  cards  adhering  to  a 
user-supplied  format. 


a  free  format  field  option  is  available  that  allows  the  user  to 
separate  input  values  by  blanks  or  commas  rather  than  using  the  standard 
(default  =  *FIXED)  10  fields  of  8  columns  each.  This  option  is  activated 
by  inserting  a  separate  card  (preceding  the  T1  card)  that  has  *FREE  in 
columns  1-5.  If  subsequent  jobs  are  provided  in  standard  format,  a  *FIXED 
card  MUST  be  inserted  before  the  T1  card  of  the  affected  data  set.  The 
♦FREE  format  option  DOES  NOT  affect  the  QQ,  H,  N,  daily  GETUSGS 
streamflow,  monthly  GETUSGS  streamflow,  or  EQ  cards.  These  cards  must  all 
be  in  the  exact  format  described  in  the  subsection  for  each  card. 


Printing  of  the  original  input  data  (including  flow  values)  is 
possible  by  use  of  the  *N0PRINT  and  *PRINT  cards,  respectively.  These  two 
cards  may  be  used  throughout  the  job  input  data  except  that  their  use  must 
not  interrupt  the  sequence  of  streamflow  data  (QQ  cards).  Whenever  one  of 
these  four  commands  (*FIXED,  *FREE,  *PRINT,  *N0PRINT)  is  provided,  it 
remains  in  effect  until  its  logical  counterpart  is  encountered.  Jobs 
following  a  seasonal  flow  duration  analysis  will  not  be  processed. 
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HYDUR  INPUT  DESCRIPTION 


The  different  values  a  variable  may  assume  and  the  conditions  for 
each  are  described  for  each  variable.  Some  variables  are  used  simply  to 
indicate  whether  or  not  a  program  option  is  to  be  used.  The  values  for 
these  variables  are  integers  and  must  be  right  justified  (punched  on  the 
far  right  side  of  the  field)  without  any  decimal  points.  Other  variables 
are  assigned  numbers  which  express  the  variables'  magnitude.  For  these, 
either  a  or  a  "+  or  is  shown  in  the  description  under  "value"  and 
the  numerical  value  of  the  variable  is  entered  as  input.  A  "+"  indicates 
only  a  positive  numerical  value  for  this  variable  is  appropriate.  The  "+ 
or  allows  for  any  non-zero  value  to  be  assigned,  where  the  variable's 
value  is  to  be  zero,  the  field  may  be  left  blank,  since  a  blank  field  is 
read  as  zero  and  any  number  without  a  sign  is  considered  positive.  Unless 
noted  othei'wise,  variable  names  beginning  with  the  letters  I,  J,  K,  L,  M, 
or  N  represent  integer  variables  and  a  decimal  point  must  not  appear  in 
the  field.  Certain  variables  (see  the  FL  card,  Field  1)  in  the  Users 
Manual  require  that  an  exact  alphanumeric  code  is  to  be  entered  to 
indicate  a  program  option.  The  values  should  be  ri jnt  justified  in  the 
indicated  field  of  the  card  without  quotations.  The  location  of  variables 
on  cards  is  sometimes  referred  to  by  an  abbreviated  designation,  for 
example,  FM.4  means  the  fourth  field  of  the  FM  card. 


All  cards  are  optional  unless  they  are  specifically  identified  as 
being  required.  The  default  values  for  each  variable,  if  assigned  by  the 
program,  are  shown  in  parentheses  under  the  variable  name. 


Several  jobs  may  be  processed  at  the  same  time  by  stacking  the 
respective  data  decks.  The  program  automatically  saves  the  values  from 
each  input  card.  Subsequent  jobs  only  require  new  cards  or  cards  whose 
values  have  changed  from  the  previous  job.  Whenever  an  input  card 
replaces  a  card  from  the  previous  job,  all  input  fields  on  the  new  card 
must  be  specified.  Caution:  If  a  subsequent  job  does  not  need  a 
particular  card  that  was  used  in  a  previous  job,  then  a  card  with  the 
proper  2-column  identifier  should  be  included  with  all  fields  left  blank. 
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HYDUR  INPUT  DESCRIPTION 


Summary  of  Program  Usage  Options* 

OPTION  DESCRIPTION  NECESSARY  CARDS 

1A  Calculation  of  a  flow  T1,  T2,  T3,  T4,  FL,  EJ,  plus  one  of 
duration  curve  from  the  following: 

streamflow  data  (1)  QQ  cards  including  sets  of  H,  N 

the  37  daily  flow  cards  for  each 
year 

(2)  or  record;  H  and  N  cards  and 
monthly  flows; 

(3)  or  flows  in  user  specified 
format  (FM  card  required). 

IB  Calculation  of  seasonal  Tl,  T2,  T3,  T4,  SD,  FL,  EJ  plus  one 
duration  curves (s)  of  the  following: 

from  streamflow  data  (1)  QQ  cards  including  sets  of  H,  N 

and  37  daily  flow  cards  for  each 
year  of  record; 

(2)  H  and  N  cards  and 
monthly  flows; 

(3)  or  flows  in  user  specified 
format  (FM  card  required). 

2A  Calculation  of  average  Cards  from  Option  1  and  at  users's 
annual  energy  from  discretion  the  PD,  PQ,  TQ,  TE,  TW, 

streamflow  data  TH  card(s). 

2B  Calculation  of  seasonal  Same  as  Option  IB  with  the  addition 
energy  production  of  SI,  S2,  SO  and/of  SS  cards, 

from  streamflow  data 

3  Calculation  of  average  Tl,  T2,  T3,  T4,  FD,  EJ,  plus  one  of 
annual  energy  from  the  following: 

flow  duration  data  (1)  EQ  cards; 

(2)  streamflow  duration  coordinates 
in  GETUSGS  format;  or 

(3)  streamflow  duration  coordinates 
in  user  supplied  format  (FM 
card  required). 

As  options  the  PD,  PQ,  TQ,  TE,  TH 

card(s ) . 


‘Required  cards  for  all  options:  Tl,  T2,  T3,  T4,  EJ  either  FL  or  FD  must 
be  provided.  PD  is  usually  a  required  card. 
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HYDUR  INPUT  DESCRIPTION 


Summary  of  Program  Usage  Options* 


OPTION  DESCRIPTION 


NECESSARY  CARDS 


4  Optimization  of 
capacity 


Option  2  or  Option  3  cards  plus  the 
OC  card.  If  costs  and  benefits  are 
analyzed  these  data  must  be  entered 
as  explained  in  Option  5  and  6. 


5 


Calculation  of  project 
benefits 


Either  the  PB  card,  or  the  CB  and 
EB  cards .  Or  MC  and  ME  cards . 


6 


Calculation  of  project  Cl  and/or  C2  card3.  CF  card  is 
costs  optional. 


♦Required  cards  for  all  options:  T1,  T2,  T3,  T4,  EJ  either  FL  or  FD  must 
be  provided.  PD  is  usually  a  required  card. 
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HYDUR  INPUT  DESCRIPTION 


TABLE  OF  CONTENTS 


GENERAL  PURPOSE 

CARDS 

1 

1.1 

T1 ,  T2,  T3,  T4 

Title  cafrds  (required) . 

.  7 

1.2 

AF  CARD 

Adjustment  factors  (optional)  . 

.  8 

1.3 

C  CARD 

Comment  cards  (otpional)  . 

.  9 

STREAMFLOW  DATA 

CARDS 

2.1 

SD  CARD 

Seasonal  duration  parameters  (optional) 

.  10 

2.2 

SI  CARD 

Seasonal  upstream  diversions  (optional) 

• 

.  13 

2.3 

S2  CARD 

Seasonal  losses  at  dam  site  (optional)  .  . 

• 

.  14 

2.4 

SO  CARD 

Seasonal  operation  of  plant  (optional)  .  . 

• 

.  14 

2.5 

SS  CARD 

Seasonal  power  storages  (optional)  .... 

• 

.  15 

2.6 

FL  CARD 

Flow  card  (required  unless  FD  is  provided) 

» 

.  16 

2.7 

FM  CARD 

User  specified  format  card  (optional)  .  . 

• 

.  18 

2.8 

QQ  CARDS 

Streamflow  data  cards  ....  . 

• 

19 

FLOW 

DURATION  DATA 

3.1 

FD  CARD 

Flow  duration  card  (required  unless  FL  is 
provided)  . 

• 

3.2 

FM  CARD 

User  specified  format  card  . 

* 

3.3 

EQ  CARD 

Exceedance  flow  data  card  . 

• 
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HYDUR  INPUT  DESCRIPTION 


TABLE  OF  CONTENTS  (continued) 


4.  POWER  DATA 

4.1  PD  CARD  Power  design  parameters  card  (usually  a 

required  card)  .  30 

4.2  PQ  CARD  Power  flow  card . 33 

4.3  TQ  CARD  Discharges  values  card . 35 

4.4  TE  CARD  Table  of  efficiencies  card . 35 

4.5  TH  CARD  Table  of  headwater  elevations  card . 36 

4.6  TW  CARD  Table  of  tailwater  elevations  card . 36 

5.  POWER  BENEFIT  INFORMATION 

5.1  PB  CARD  Power  benefit  card . 37 

5.2  CB  CARD  Capacity  benefit  cards . 40 

5.3  EB  CARD  Energy  benefit  cards . 41 

5.4  MC  CARD  Seasonal  capacity  benefits . 42 

5.5  ME  CARD  Seasonal  energy  benefits . 42 
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HYDUR  INPUT  DESCRIPTION 


1  GENERAL  PURPOSE  CARDS 


T1 

T2 

T3 

T4 


1.1  T1,T2,T3  AND  T4  -  Title  Cards  (Required) 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

T1,T2 

Card  identification  characters. 

T3  or  T4 

1-10 

TITLE 

ALPHA 

Alphanumeric  information  to  identify  the 

job.  Information  on  the  T4  card  is  not 
retained  for  future  jobs. 


EXHIBIT  7 


Page  7  of  54 


AF 


HYDUR  INPUT  DESCRIPTION 


1.2  AF  CARD  -  Adjustment  Factors  {Optional) 


The  AF  card  specifies  units  of  measurement  used  in  HYDUR. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

AF 

Card  identification  characters. 

1 

SYSTEM 

(ENGLSH) 

METRIC 

All  input  values  are  in  metric  units. 

Flows  in  cubic  meters/second  and 
lengths  in  meters . 

ENGLSH 

All  input  values  are  in  english  units. 

Flows  in  cubic  feet/second  (cfs)  and  lengths 
in  feet. 

2 

ALEN 

(1.0) 

+ 

Adjustment  factor  multiplied  times  the 
display  of  all  length  (feet  or  meters) 
related  parameters. 

3 

APWR 

(1.0) 

+ 

Adjustment  factor  multiplied  times  the 
display  of  all  power  (kw)  related 
parameters . 

4 

AENRG 

(1.0) 

+ 

Adjustment  factor  multiplied  times  the 
display  of  all  energy  (mwh)  related 
parameters . 

5 

ADOLR 

(1.0) 

+ 

Adjustment  factor  multiplied  times  the 
display  of  all  money  ($)  parameters. 

6 

AFLOW 

(1.0) 

+ 

Adjustment  factor  multiplied  times  the 
display  of  all  flow  (cfs  or  cats)  related 

parameters . 
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1.3  C  CARDS  -  Comment  Cards  (Optional) 


Comment  cards  enable  the  user  to  insert  documentation  information 
throughout  the  input  card  stream.  Comment  cards  are  displayed  in  the  echo 
printing  of  the  input  job  (T1-EJ)  but  have  no  effect  on  the  operation  of 
the  program.  Any  number  of  comment  cards  can  be  inserted.  The  only 
restriction  is  that  they  cannot  be  inserted  between  flow  values  (FL 
through  QQ  cards ) . 
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HYDUR  INPUT  DESCRIPTION 


2  STREAMFLOW  DATA  CARDS  (Optional) 


Streamflow  cards  SD,  SI,  S2,  SO,  SS,  FL,  FM,  QQ,  are  used  to  develop 
a  duration  curve  from  streamflow  values.  The  streamflow  cards  in  this 
section  are  omitted  if  the  streamflow  duration  curve  is  input  directly  or 
if  a  flow  duration  curves  is  used  from  a  previous  job  in  stacked  jobs. 
These  cards  are  required  only  when  flow  duration  cards  are  omitted. 


2.1  SD  CARD  -  Seasonal  Duration  Parameters  (Optional) 

The  user  is  cautioned  that  the  seasonal  duration  option  of  HYDUR  has 
not  undergone  extensive  testing.  If  probelms  occur,  it  is  recommended 
that  the  normal  annual  analyses  be  performed. 


The  SD  card  is  used  to  develop  streamflow  duration  curves  for 
specified  time  intervals.  The  time  interval  (season)  may  vary  from  one 
month  to  one  year  in  duration.  This  card  is  used  only  in  conjunction  with 
streamflow  values  input  on  FL  cards.  The  SD  card  must  precede  the  FL  card 
in  the  input  sequence  of  a  job.  If  nc  SD  card  is  supplied  or  if  the  above 
requirement  is  not  met,  an  annual  streamflow  duration  curve  is  developed. 


FIELD 

VARIABLE  VALUE 

DESCRIPTION 

0 

ID  SD 

Card  identification  characters. 

1 

INTBIAS  + 

(0) 

The  number  of  streamflow  values  to  skip 
before  beginning  the  first 

streamf low-duration  analysis.  Maximum  value 
is  defined  in  Field  3  (NINYR). 

2 

NOLUDE  + 

(value  in  SD.3) 

Defines  the  number  of  consecutive 
streamflow  values  to  include  in  the 

development  of  duration  curve (s)  once  the 
proper  starting  position  has  been 
determined.  Additional  flows  from 
subsequent  years  of  record,  which  are  to  be 
incorporated  in  the  duration  curve,  are 
found  by  skipping  ((SD.3)  -  (SD.2))  values. 
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FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

3 

NINYR 

(0) 

+ 

A  value  is  required  if  seasonal  analysis  is 
to  be  performed.  Defines  the  number  of 
streamflow  values  which  comprise  a  year 

For  monthly  , streamflow  records,  this  value 
would  be  twelve  (12).  For  daily  flow 
records,  this  value  should  be  365  for 
blocked  records  required  (12*31)  or  372 
daily  flows  to  be  read  per  year.  This  type 
of  input  requires  negative  flow  values  to  be 
read  for  non-existent  days,  such  as  February 
31.  If  FL.2  =  DAILY  and  NINYR=372,  then 
GETUSGS  daily  blocked  data  will  be  expected. 

4  NSEASN  +  Defines  the  number  of  seasons  to  be 

(1)  analyzed.  A  different  streamflow  duration 

curve  is  developed  for  each  season. 

The  number  of  streamflow  values  to  skip  for 
each  season.  This  value  is  added  to  an 
accumulator  each  time  a  new  season  is 
processed  and  the  accumulated  value 
determines  the  starting  position  for  each 
season.  The  user  can  set  this  parameter  to 
allow  for  overlapping  of  streamflows  from 
one  season  to  the  next.  The  maximvuu  value 
for  this  parameter  is  determined  by  the 
ratio  ( (SD.3)/(SD.4)) . 

or  0  For  annual  power  analyses  (SD.2  =  SD.3),  no 
restriction  to  an  upper  limit  is  considered. 
For  seasonal  power  analyses  (SD.2  less  than 
SD.3),  the  program  selects  this  upper  limit 
based  on  scanning  the  flow  record  (QQ 
cards ) . 


*The  use  of  this  option  results  in  the  OP  card  being  ignored.  Instead, 
intermediate  capacities  are  determined  based  U£-on  a  constant  interval 
determined  by  the  following  equation: 

BINTV  =  ( C AP YMX-C APYMN )/ (N— 1 ) 

where : 

N  represents  the  total  number  of 
capacity  values  to  consider. 


6*  CAPYMX 

(0) 


5  NSEASKP  l- 

( value  in  SD.2) 
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HYDUR  INPUT  DESCRIPTION 


FIELD  VARIABLE  VALUE 


DESCRIPTION 


7*  CAPYMN 
(0) 


+  The  maximum  installed  capacity  (in 

kilowatts)  considered  in  subsequent 

optimization  analyses . 

-  or  0  For  annual  power  analyses  (SD.2  =  SD.3),  no 
restriction  to  a  lower  limit  of  installed 
capacity  in  subsequent  optimization  analyses 
is  donsidered.  For  seasonal  power  analyses 
(SD.2  less  than  SD.3),  the  program  selects 
this  lower  limit  based  on  scanning  the  flow 
record  (QQ  cards). 

+  The  minimum  installed  capacity  (in 

kilowatts)  considered  in  subsequent 

optimizations  analyses. 


♦The  use  of  this  option  results  in  the  OP  card  being  ignored.  Instead, 
intermediate  capacities  are  determined  based  upon  a  constant  interval 
determined  by  the  following  equation: 

BIN TV  =  ( CAPYMX-CAPYMN )/ (N-1 ) 

where : 

N  represents  the  total  number  of 
capacity  values  to  consider. 
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The  following  set  of  optional  cards  (SI,  S2,  SO  and  SS)  enable  the 
user  to  vary  parameters  which  typically  change  with  time  (season).  An  SD 
card  must  be  input  in  this  job. 


2.2  SI  CARDS  -  Seasonal  Upstream  Diversions  (Optional) 

The  Si  card  specifies  upstream  diversions  to  be  varied  by  season. 
The  use  of  this  card  results  in  variable  DIV  (PQ.3)  to  be  overridden  for 
each  season  analyzed.  A  maximum  of  12  values  may  be  input  (two  SI  cards); 
however,  only  the  first  SD.4  values  will  be  used  in  the  concurrent  job. 


FIELD  VARIABLE  VALUE  DESCRIPTION 


0  ID  SI  Card  identification  characters. 

1-10  SDIV  +  or  -  Upstream  diversions  (in  cfs  or  cms)  varying 

by  season.  Up  to  12  values,  input  in  groups 
of  10  values  per  card,  are  allowed.  Once 
the  maximum  number  has  been  defined, 
additional  SI  card  input  is  ignored. 
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HYDUR  INPUT  DESCRIPTION 


2.3  S2  CARDS  -  Seasonal  Losses  at  Dam  Site  (Optional) 


The  S2  card  specifies  seasonable  tailwater  losses.  The  S2  card 
overrides  variable  QLOSS  (PQ.4).  Input  limitations  and  requirements  are 
identical  to  those  defined  for  the  SI  card. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

S2 

Card  identification  characters. 

1-10 

SQL 

+  or  - 

Seasonal  losses  effecting  tailwater  (in 
cfs  or  cms). 

2.4  SO  CARDS  -  Seasonal  Operation  of  Plant  (Optional) 


The  SO  card  specifies  seasonal  plant  operation  (seasonal  plant 
factor).  The  SO  card  overrides  variable  VAPP  (PD. 6).  Input  limitations 
and  requirements  are  identical  to  those  defined  for  SI  card. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

SO 

Card  identification  characters. 

1-10 

SPF 

+ 

Plant  factor  (greater  than  zero  and  less 
than  one)  values. 
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HYDUR  INPUT  DESCRIPTION 


ss 


2.5  SS  CARDS  -  Seasonal  Power  Storages  (Optional) 


The  SS  card  specifies  seasonal  reservoir  storage  (power  storage 
ratios).  The  SS  card  overrides  variable  PSR  (PD. 7).  Input  limitations 
and  requirements  are  identical  to  those  defined  for  the  SI  card. 

FIELD  VARIABLE  VALUE  DESCRIPTION 


0  ID  SS  Card  identification  characters. 

1-10  SSR  Seasonal  power  storage  ratios  (greater  than 

or  equal  to  zero). 
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HYDUR  INPUT  DESCRIPTION 


2.6  FL  CARD  -  Flow  Card  (Required  unless  FD  is  provided) 


The  FL  card  is  required  if  streamflow  data  are  used  to  calculate  the 
streamf low  duration  curve.  This  card  is  omitted  when  a  streamflow 
duration  curve  is  input  on  FD  cards  or  when  a  flow  duration  curve  from  a 
previous  job  is  used. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

FL 

Card  identification  characters. 

1 

FTYPE 

(card) 

CARD 

Streamflow  data  will  be  read  by  cards 
(input  file). 

FILE 

Streamflow  data  will  be  read  from  a  separate 
alternative  file  on  Unit  3.  (See  Exhibit  6 
in  Users  Manual . ) . 

2 

FMTYPE 

(STD) 

STD 

Streamflow  values  are  provided  on  QQ  cards 
in  standard  card  format  of  10F8.0. 

DAILY 

Average  daily  streamflow  data  will  be 
provided  in  GETUSGS  sequential  or  blocked 
format  including  header  cards. 

MONTH 

Average  monthly  streamflow  data  will  be 
provided  in  GETUSGS  format  including  header 
cards . 

USER 

Streamflow  data  is  provided  in  user  supplied 
format.  When  this  option  is  selected,  the 
FM  card  must  be  the  next  input  card 
provided. 

3 

PUNCH 

(NOPUNCH) 

NOPUNCH 

Streamflow  duration  ordinates  will  not  be 
punched  on  cards  during  this  job. 

PUNCH 

Pairs  of  percent  of  time  and  streamflow 
ordinates  on  the  streamflow  duration  curve 
will  be  punched  in  standard  (FL.2  equals 
STD)  format  of  10F8.0  during  this  job. 

4 

POWANAL 

(YES) 

YES 

Power  analysis  will  be  performed  (after 
generating  the  streamflow  duration  curve). 

PD  card  is  required. 
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HYDUR  INPUT  DESCRIPTION 


FIELD 


5 


VARIABLE  VALUE 


DESCRIPTION 


LENREC 

(80) 


FLOWLO 

(0) 


NO  Power  analysis  will  not  be  performed.  Only 

streamf low-duration  curve  will  be  generated. 

This  parameter  is  required  if  variable 
FTYPE  (FL.1)  equals  FILE.  FILE  defines  the 
record  length  of  the  alternative  file. 

0  Record  length  is  assumed  to  be  equal  to  a 

typical  card  image  of  eighty  (80)  columns. 

+  Record  length  of  alternate  file  (not  to 

exceed  150).  For  example,  if  streamf lows  on 
the  alternative  file  are  written  using  the 
format  of  (7X,5F8.2),  then  the  corresponding 
record  length  would  be  ((7+ (5*8))  or  47. 

+,-,0  Discharge  values  less  than  this 

value  (FLOWLO)  will  not  be  considered 
in  determining  the  flow-duration  curve. 
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HYDUR  INPUT  DESCRIPTION 


2.7  FM  CARD  -  User  Specified  Format  Card  (Optional) 

The  FM  card  is  required  when  variable  FMTYPE  (FL.2  or  FD.2)  is  equal 
to  USER,  the  FM  card  indicates  format  of  the  user  supplied  data  and  the 
number  of  streamflow  or  streamflow  duration  data  values.  Either  variable 
TERM  ( FM . 1 )  which  specifies  a  user  supplied  terminator  code,  or  variable 
ITVALS  ( FM . 2 )  which  specifies  the  number  of  streamflow  data  values  or 
pairs  of  streamflow  duration  ordinajtes,  must  be  provided. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

FM 

Card  identification  characters. 

1 

TERM 

(-999) 

0 

Variable  ITVALS  (FM.2)  will  be  used  to 
indicate  the  number  of  streamflow  values 
or  pairs  of  streamflow  duration  data  values 
or  the  default  value  (-999)  will  be  used  if 
ITVALS  is  blank. 

+  or  - 

Streamflow  or  streamflow  duration  ordinates 
are  to  be  read  until  TERM  is  encountered. 
The  absolute  value  of  TERM  should  be  greater 
than  or  equal  to  1.0,  and  placed  in  the 
field  immediately  following  the  last  data 
value . 

2 

ITVALS 

(0) 

0 

Variable  TERM  (FM.1)  indicates  the  end  of 
the  streamflow  or  flow  duration  data. 

+ 

Number  of  user  specified  streamflows  (or 
pairs  of  flows  and  percent  of  time  exceeded 
ordinates)  to  be  read. 

3 

NVALS 

+ 

Number  of  streamflow  (or  pair)  values 
specified  in  the  user  designed  format  on 
fields  4-10  of  the  FM  card. 

4-10 

FMT 

ALPHA 

User  supplied  format  for  streamflow  data  or 

the  flow  duration  curve.  Flow  duration 
ordinates  must  be  specified  in  pairs  of 
percent  of  time  exceeded  and  flow  values . 
Parentheses  must  be  included,  e.g., 

( 8X , 6F6 . 0 ) .  The  specification  may  be 
continued  to  the  following  card  by 
specifying  another  FM  card  and  continuing 
the  format  in  column  3. 
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QQ 


2.8  QQ  CARDS  -  Streamflow  Data  Cards 


The  QQ  cards  specify  streamflow  values  when  variable  FMTYPE  (FL.2)  is 
equal  to  "STD".  The  QQ  is  not  required  in  the  first  two  columns.  The 
last  streamflow  value  must  be  -999.  The  QQ  cards  follow  the  FL  card. 


FIELD  VARIABLE  VALUE  DESCRIPTION 


1-10  FLOWS  +  Streamflow  values  in  cfs  or  cms . 

Missing  data  values  should  be  indicated  by 
entering  a  negative  number  (see  Field  FL.6). 

-999  Required  terminator  code  that  must 

immediately  follow  the  last  streamflow 
value . 
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FD 


3  FLOW  DURATION  DATA  (Required  unless  streamflow  data  cards  are  used) 


The  FD,  FM,  and  EQ  cards  specify  criteria  for  inputting  streamflow 
duration  curves.  When  individual  streamflow  values  are  provided  the  cards 
in  this  section  are  omitted. 


3.1  FD  CARD  -  Flow  Duration  Card  (Required  unless  FL  is  provided) 


The  FD  card  is  required  if  streamflow  duration  data  are  to  be  input. 
This  card  is  omitted  when  streamflow  data  are  specify  on  the  FL  card  or  a 
flow-duration  curve  from  a  previous  job  is  used. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

FD 

Card  identification  characters. 

1 

FTYPE 

(CARD) 

CARD 

The  flow  duration  curve  will  be  input 
cards . 

by 

FILE 

The  flow  duration  curve  will  be  input 
separate  alternative  file  on  Unit  3. 
Exhibit  6  in  Users  Manual) . 

from  a 
(See 

2 

FMTYPE 

STD 

Streamflow  duration  and  corresponding 

percent  of  time  exceeded  values  are  provided 
in  standard  card  format  -  10F8.0  (see  EQ 
cards).  Data  must  be  terminated  by  -999. 


USER  Streamflow  duration  and  percent  of  time 

exceeded  values  are  provided  in  user 
supplied  format.  The  FM  card  must 
immediately  follow  the  FD  card. 
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FD 


TABLE  2:  Standard  Two  Letter  Postal  and  Two  Digit  FIPS  Codes 


AL 

— 

01 

— 

ALABAMA 

AK 

— 

02 

— 

ALASKA 

AZ 

— 

04 

— 

ARIZONA 

AR 

— 

05 

— 

ARKANSAS 

CA 

— 

06 

— 

CALIFORNIA 

CO 

— 

08 

— 

COLORADO 

CT 

— 

09 

— 

CONNECTICUT 

DE 

— 

10 

— 

DELAWARE 

DC 

— 

11 

— 

DISTRICT  OF  COLUMBIA 

FL 

— 

12 

— 

FLORIDA 

GA 

— 

13 

— 

GEORGIA 

HI 

— 

15 

— 

HAWAII 

ID 

— 

16 

— 

IDAHO 

IL 

— 

17 

— 

ILLINOIS 

IN 

— 

18 

— 

INDIANA 

IA 

— 

19 

— 

IOWA 

KS 

— 

20 

— 

KANSAS 

KY 

— 

21 

— 

KENTUCKY 

LA 

— 

22 

— 

LOUISIANA 

ME 

— 

23 

— 

MAINE 

MD 

— 

24 

— 

MARYLAND 

MA 

— 

25 

— 

MASSACHUSETTS 

MI 

— 

26 

— 

MICHIGAN 

MN 

— 

27 

— 

MINNESOTA 

MS 

— 

28 

— 

MISSISSIPPI 

MO 

— 

29 

— 

MISSOURI 

MT 

— 

30 

— 

MONTANA 

NE 

— 

31 

— 

NEBRASKA 

NV 

— 

32 

— 

NEVADA 

NH 

— 

33 

— 

NEW  HAMPSHIRE 

NJ 

— 

34 

— 

NEW  JERSEY 

NM 

— 

35 

— 

NEW  MEXICO 

NY 

— 

36 

— 

NEW  YORK 

NC 

— 

37 

— 

NORTH  CAROLINA 

ND 

— 

38 

— 

NORTH  DAKOTA 

OH 

— 

39 

— 

OHIO 

OK 

— 

40 

— 

OKLAHOMA 

OR 

— 

41 

— 

OREGON 

PA 

— 

42 

— 

PENNSYLVANIA 

RI 

— 

44 

— 

RHODE  ISLAND 

SC 

— 

45 

— 

SOUTH  CAROLINA 

SD 

— 

46 

— 

SOUTH  DAKOTA 

TN 

— 

47 

— 

TENNESSEE 

TX 

— 

48 

— 

TEXAS 

UT 

— 

49 

— 

UTAH 

VT 

— 

50 

— 

VERMONT 

VA 

— 

51 

— 

VIRGINIA 

WA 

— 

53 

— 

WASHINGTON 

WV 

— 

54 

— 

WEST  VIRGINIA 

WI 

— 

55 

— 

WISCONSIN 

WY 

— 

56 

— 

WYOMING 
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3.2  FM  CARD  -  User  Specified  Format  Card 


The  FM  card  is  required  to  specify  the  format  of  the  duration  data 
when  FD.2  equals  "USERS"  (see  Section  2.7). 


3.3  EQ  CARD  -  Exceedance  Flow  Data  Card 


The  EQ  card  is  required  to  specify  ordinates  on  the  streamflow 
duration  curve  when  FD.2  is  equal  to  "STD".  The  "EQ"  is  not  required  in 
the  first  two  columns.  The  last  pair  of  streamflow  duration  coordinates 
must  be  followed  by  a  -999.  A  maximum  of  70  pairs  of  PE  and  QQ  values  may 
be  input. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

EQ 

Card  identification  characters  (optional). 

1 

PE  ( 1 ) 

+ 

First  ordinate  for  percent  of  time  exceeded 
oridnate  expressed  as  percent  or  as  a 
decimal  fraction.  Values  input  in 

decreasing  order. 

-999 

Required  terminator  code  that  immediately 
follows  the  last  streamflow  value. 

2 

QQ(1) 

+ 

First  ordinate  streamflow  in  cfs  or  cubic 
meters/second  corresponding  with  variable  PE 
(EQ.1),  values  are  input  in  increasing 
order . 

3-10 

PE(N)fQQ(N) 

+ 

Same  as  above  for  up  to  70  pairs  of 

ordinates . 
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GETUSGS 


3.4  GETUSGS  STREAMFLOW  DURATION  CURVE 


These  cards  are  required  when  FD.2  is  equal  to  "FLDUR" .  The  GETUSGS 
consists  of  a  header  card  followed  by  pairs  of  flow  duration  ordinates. 
The  number  of  pairs  of  streamflow  duration  coordinates  must  be  indicated 
by  either  "NPQ"  on  the  header  card  or  the  last  percent  of  time  ordinate 
must  be  followed  by  a  -999.  Must  follow  the  FD  card. 


NOTE:  GETUSGS  data  is  stored  in  english  units  only  (cfs).  Therefore, 
using  the  metric  option  requires  an  adjustment  factor  to  be  specified  in 
field  (PQ.5 ) . 


HEADER 

CARD 

FORMAT 

COLUMN 

VARIABLE 

VALUE 

DESCRIPTION 

12 

1-2 

ISTATE 

+ 

Two  digit  state  FIPS  code  from 
Table  2  (Exhibit  7). 

18 

3-10 

IGAUGE 

+ 

Eight  digit  integer  USERS  gage 
number . 

12 

11-12 

NPQ 

+ 

Number  of  pairs  of  discharge  and 
Percent  of  time  exceeded  ordinates 
used  to  describe  the  streamflow 
duration  curve. 

DATA  CARDS 


Supply  data  cards  as  needed  -  5  pairs  of  (QQ,  PE)  per  card. 


FORMAT 

COLUMN 

VARIABLE 

VALUE 

DESCRIPTION 

F10.2 

1-10 

QQ 

+ 

First  discharge  in  cfs  (lowest 
discharge) . 

F6.1 

11-16 

PE 

+ 

First  percent  of  time  exceeded 
ordinate  corresponding  to  the 
discharge. 

4(F10.2,F6. 1 ) 

17-80 

QQ/PE 

+ 

Subsequent  pairs  of  QQ  and  PE 

values  are  provided. 
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4  POWER  DATA 


The  PD,  PQ,  TQ,  TE,  TH,  and  TW  cards  define  the  power  generating 
characteristics  of  the  power  plant,  the  assumed  head,  and  flow 
characteristics , 


4.1  PD  CARD  -  Power  Design  Parameters  Card  (usually  a  required  card) 

The  PD  card  specifies  design  parameters  for  the  power  facility.  The 
value  CAPDES  specifies  if  an  optimization  analysis  is  to  be  performed. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

PD 

Card  identification  characters. 

1 

QDES 

0 

The  maximum  penstock  discharge  capacity  will 

be  calculated  based  on  the  installed 
capacity,  overload  factor  of  the  plant, 
efficiency,  and  head  (see  Equation  2).  For 
optimization  analyses,  variable  CAPDES 
(PD. 4)  equals  zero,  QDES  is  set  equal  to  the 
selcted  discharge  from  the  flow  duration 
curve . 

+  Maximum  penstock  discharge  capacity  in  cfs 

or  cms . 

2  EFF  0  Powerplant  efficiency  (TE  card) 

(0.86)  corresponding  to  the  maximum  penstock 

discharge  capacity  will  be  used.  If  TE  card 
is  not  input  the  efficiency  of  .86  is  used. 

+  The  design  powerplant  efficiency  expressed 

as  a  decimal  fraction.  Powerplant 
efficiency  equals  the  efficiency  of  the 
turbine  times  the  efficiency  of  the 
generator . 

3  HEAD  0  The  head  will  be  determined  from  the  head 

elevation  (TH  card)  and  tailwater  elevation 
(TW  card)  data. 
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FIELD 


4 


5 


6 


VARIABLE  VALUE 


DESCRIPTION 


+ 


CAPDES  0 

(0) 


+ 


OVLOAD  + 

(1.15) 


UAPF 

(1.0) 


0 


+ 


PSR  0 


The  net  design  head  of  the  turbine  in  feet 
or  meters.  Constant  head  for  all  power 
calculations,  head  elevation  tailwater  data 
(TH)  and  TW  are  not  input. 

Plant  capacity  will  be  optimized  based  on 
criteria  specified  on  the  OC  card. 

Installed  powerplant  nameplate  capacity  in 
kilowatts . 

Overload  factor  for  the  power  installation. 
Input  as  a  decimal  to  be  multplied  by 
variable  CAPDES  (PD. 4). 

User  defined  annual  plant  factor  which 
represents  the  annual  plant  factor  required 
to  meet  annual  firm  energy  requirements. 
UAPF  can  be  expressed  as: 

UAPF  =  AFE/(DC*8.76) 

•  ere: 

AFE  =  annual  firm  energy  in  mwh 
DC  =  dependable  capacity  in  kw 

A  value  other  than  1.0  will  affect  net  power 
head  calculations  since  the  factor  will  be 
used  to  adjust  discharge  values  to 
approximate  the  flows  at  the  time  of  actual 
hydropower  operation. 

Annual  plant  factor  value  is  not  used  to 
adjust  flows  in  the  calculation  of  the 
tailwater  elevation. 

Annual  plant  factor  (value  between  0  and  1), 
divided  into  the  discharge  values  to  account 
for  increased  tailwater  elevations 

encountered  in  peaking  operations. 

No  adjustment  to  flow-duration  curve.  Flow 
duration  curve  assumed  to  represent  flow 
pattern  available  for  power  generation. 
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HYDUR  INPUT  DESCRIPTION 


FIELD 


8 


VARIABLE  VALUE 


DESCRIPTION 


+ 


XUNITS  0 


-t- 


PEAKF  0 

(1.0) 


+ 


Power  storage  ratio  (PSR)  used  to  adjust  the 
flow-duration  curve  to  account  for  project 
storage  (see  Exhibit  2).  Equals  the  power 
storage  in  acre-feet  divided  by  the  product 
of  722.50  and  average  annual  inflow  in  cfs. 

Power  storage  ratio  (PSR)  used  to  adjust  the 
flow-duration  curve  to  account  for  project 
storage  as  a  function  of  the  average  annual 
draft  rate  (see  Exhibit  2). 

The  number  of  units  required  to  deliver  the 
plant’s  installed  capacity  will  be 
determined  by  selecting  a  unit  design 
discharge  equal  to  the  flow  corresponding  to 
an  exceedence  percent  of  75. 

EXAMPLE:  The  design  discharge  of  a  plant  is 
5000  cfs  and  that  the  75  percent  exceedence 
flow  value  is  3000  cfs,  the  number  of  units 
assigned  by  the  program  would  be: 

( 5000/3000 )+1  =  2.67  or  2  units 

Number  of  units  (up  to  50)  assumed  in  place 
for  this  job.  Combined  capacities  equal  the 
total  installed  plant  capacity. 

NOTE:  If  the  project  is  presently  in  place, 
and  ECAP  (Cl. 4)  is  greater  than  zero, 
existing  and  total  units  assumed  installed 
are  defined  as:  existing  units  to  the  right 
of  the  decimal  and  total  units  to  the  left . 
For  example,  the  input  value  to  define  a 
total  of  four  (4)  units,  three  (3)  of  which 
are  presently  in  place,  would  be  4.03. 

Peaking  adjustin'  c.  factor.  No  adjustment 
to  capacities  determined  in  the  optimization 
analysis  will  be  made. 

The  peaking  adjustment  factor  (between  1  and 
10)  will  be  used  to  adjust  (multiplied  by) 
design  discharge  values  obtained  in  the 
optimization  analysis.  NOTE:  If  the  maximum 
installed  capacity,  variable  CAPYMX  (SD.6), 
is  greater  than  zero  no  adjustment  is  made. 
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4.2  PQ  CARD  -  Power  Flow  Card 


The  PQ  card  specifies  flow  constraints 
produce  power. 


on  the  water  available  to 


FIELD  VARIABLE  VALUE  DESCRIPTION 


0  ID 

1  QSUB 

(0) 


2  QMIN 

(0) 


DIV 

(0) 


PQ  Card  identification  characters. 

0  Submergence  flow  is  not  considered. 


+ 


+ 


Submergence  flow  in  cf s  or  cms .  No  power 
generated  for  discharge  values  greater  than 
input  value. 


Minimum  flow  in  cfs  or  cms.  No  power  is 
generated  for  discharges  less  than  QMIN. 
Absolute  input  values  less  than  one  will 
cause  QMIN  to  be  defined  as  a  percentage  of 
the  design  discharge,  using  the  following 
equation: 


QMIN  =  PER*QD/XUNITS 

where: 

PER  =  input  value  where  absolute 

value  is  less  than  one 

QD  =  current  design  q  value. 

(this  value  will  vary  during 
optimization  runs 

XUNITS  =  Number  of  power  units 

installed  (PD. 8) 


NOTE:  QMIN  is  adjusted  to  reflect  peaking 
operation  if  variable  VAPF  (PD. 6)  is  less 
than  1.0. 


Flow  diversion  in  cfs  or  cms  above  the 
powerhouse  (average  evaporation  losses  may 
be  included  in  DIV ) .  Value  is  subtracted 
from  the  streamflow  values  on  the  flow 
duration  curve  before  performing  the  power 
analysis . 
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FIELD 

4 

5 

6 

7 


VARIABLE  VALUE 


DESCRIPTION 


QLOSST 

(0) 


QFACT 

(0) 


HMIN 

(0) 


HMAX 

(0) 


SPLEF 

(YES) 


+  Diversion  of  water  in  cfs  or  cms  around 

the  powerhouse  (fish  ladder,  leakage,  etc.). 
This  flow  is  not  used  for  power  production 
but  affects  the  headwater  and  tailwater 
elevations. 

0  No  adjustment  will  be  made  to  streamflow 

values  on  flow  duration  curve  prior  to  power 
analysis . 

+  Adjustment  (multiplier)  factor  to  streamflow 

values  on  flow  duration  curve  prior  to 
performing  the  power  analysis  (e.g., 
drainage  area  adjustment,  etc.). 

+  Minimum  head  required  to  produce  power. 

No  power  is  generated  for  powerheads 
less  th^n  input  value  (a  value  less  than  1.0 
implies  the  minimum  head  will  be  calculated 
as  the  product  of  HMIN  times  the  design 
head) . 

+  Maximum  head  allowed  to  produce  power.  No 

power  is  generated  for  powerheads  greater 
than  HMAX  (a  value  less  than  2.0  implies  the 
maximum  head  will  be  calculated  as  the 
product  of  HMAX  times  the  design  head). 

BLANK  Spillage  from  the  reservoir  will  affect 

OR  the  tailwater  conditions  at  the  powerhouse. 

YES  Typically  occurs  at  overfall  dams  where  the 
draft  tubes  exit  below  the  base  of  the  dam. 

NO  Spillage  from  the  reservoir  will  not  affect 

the  tailwater  conditions  at  the  powerhouse. 
This  condition  exists  at  installations  where 
the  powerhouse  and  spillway  are  separated 
from  one  another. 
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4.3  TQ  CARD  -  Discharges  Values  Card 


The  TQ  card  specifies  discharge  values  associated  with  the  operation 
efficiencies,  headwater  and  tailwater  elevations  (TE,  TH,  and  TW  cards, 
respectively).  A  maximum  of  10  entries  is  allowed.  TQ-TW  cards  may  be 
omitted  if  the  constant  head,  HEAD  (PD. 3)  and  efficiency,  variable  EFF 
(PD. 2)  are  input. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

TQ 

Card  identification  characters. 

1 

TQ(  1 ) 

+ 

Reservoir  discharge  in  cf s  or  cms . 
Corresponds  to  operating  efficiency  (TE(1)), 
headwater  elevation  (TH ( 1 ) ) ,  and  tailwater 
elevations  (TW ( 1 ) ) . 

2-10 

TQ(N) 

+ 

Same  as  above  for  up  to  10  values. 

4.4  TE  CARD  -  Table  of  Efficiencies  Card 


The  TE  card  specifies  the  combined  operating  efficiencies  of  the 
turbine  and  generator  units .  These  values  correspond  to  the  discharge 
values  on  the  same  field  of  the  TQ  card.  A  TQ  card  must  be  input  for  each 
TE  card. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

TE 

Card  identification  characters. 

1 

TE  ( 1  ) 

+ 

Operating  efficiencies  expressed  as  decimal 
fractions  corresponding  to  discharge  values 
on  the  first  field  of  the  TQ  card. 

2-10 

TE  (N ) 

+ 

Same  as  above  for  up  to  10  values. 

EXHIBIT  7 


Page  29  of  54 


TH 

TW 


HYDUR  INPUT  DESCRIPTION 


4.5  TH  CARD  -  Table  of  Headwater  Elevations  Card 


The  TH  card  specifies  the  elevation  of  the  headwater.  Elevations  are 
in  increasing  order.  If  this  card  is  omitted  then  the  design  head  on  the 
variable  HEAD  (PD. 3)  will  be  used  in  the  power  equation. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

IL 

TH 

Card  identification  characters. 

1 

TH  ( 1 ) 

+ 

Headwater  elevation  in  feet  or  meters 
corresponding  to  the  discharge  values  on  the 
first  field  of  the  TQ  card 

2-10 

TH  (N ) 

+ 

Same  as  above  for  up  to  10  values. 

4.6  TW  CARD  -  Table  of  Tailwater  Elevations  Card 


The  TW  card  specifies  the  elevation  of  the  tailwater. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

TW 

Card  identification  characters. 

* 

1 

TW  ( 1 ) 

+ 

Tailwater  elevatic  .  in  feet  or  meters 
corresponding  to  the  discharge  vaue  on  the 
first  field  of  the  TQ  card. 

2-10 

TW  (N ) 

+ 

Same  as  above  for  up  to  10  values. 
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5  POWER  BENEFIT  INFORMATION 


Benefit  information  cards  (PB,  CB,  EB,  MC,  ME,  and  EF  cards)  are 
required  when  benefit  computations  are  to  be  performed.  If  the 
optimization  analysis  (OC  card)  specif ieds  benefits,  net  cost,  or  B/C 
ratios,  etc.  are  to  be  optimized  then  either  variable  IREG  (PB.1),  or  (CB 
and  EB),  or  (MC  and  ME)  card  are  required. 


5.1  PB  CARD  -  Power  Benefit  Card 


The  PB  card  is  required  when  power  benefits  from  one  of  the  32  FERC 
regions  (see  Table  IB)  are  desired. 

FIELD  VARIABLE  VALUE  DESCRIPTION 


0  ID  PB 

1  IREG  0 

(0) 


+ 


2  CBR 

(0.50) 


Card  identification  characters. 

Power  and  energy  benefits  are  specified 
on  CB  and  EB  cards  or  MC  and  ME  cards. 

Two  digit  integer  region  code  from  Table  1A/B 
(page  14)  main  test.  FERC  energy  and 
capacity  benefits  from  the  region 
illustrated  on  Figure  4  in  the  Users  Manual 
are  applied  in  the  program. 

NOTE:  Benefit  relationship  developed  are 
functionally  related  to  plant  factor,  since 
this  parameter  is  often  a  good  indicator  of 
the  least-cost  alternative  to  hydropower 
facilities.  Storage  projects  often  use  this 
approach  in  determining  benefits.  However, 
it  may  not  be  appropriate  for  run-of-river 
projects.  Values  for  energy  (capacity  if 
appropriate)  should  be  determined  by  fuel 
displacement  (or  similar)  analysis  and  input 
using  the  EB  card. 

Capacity  benefit  reduction  factor  that  is 
multiplied  times  the  capacity  benefit  to 
determine  the  interruptible  capacity 
benefit. 
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FIELD 


3 


4 


5 


6 


VARIABLE  VALUE 
0 


+ 


DCAP  0 

(value  in  PB.6) 


+ 


AAEST  0 


+ 


RAAE  0 


+ 


SRP  0 

(0.85) 


+ 


ERATIO  0 


+ 


DESCRIPTION 

A  capacity  reudction  factor  of  .50  will 
be  multiplied  t  ones  the  capacity  benefit  to 
determine  the  interruptible  capacity 
benefit. 

Capacity  benefit  reduction  factor  (between  0 

and  1.0)  to  be  multiplied  times  the  capacity 

to  determine  the  interruptible  benefit. 

* 

The  dependable  capacity  is  based  on 
the  value  input  on  variable  SRP(PB.6). 

Dependable  capacity  in  kilowatts  used  to 
calculate  the  capacity  benefit. 

Average  annual  energy  will  be  calculated  by 
the  program  to  determine  the  energy  benefit. 

Average  annual  energy  in  kilowatt-hours  to 
be  used  in  calculation  of  the  energy  benefit 
(in  place  of  calculated  energy). 

No  ratio  is  used.  No  adjustment  will  be 
made  to  the  calculated  energy  used  to 
determine  benefits. 

Ratio  (decimal  vaue)  multiplied  times  the 
average  annual  energy  to  adjust  the  energy 
used  in  the  benefit  computations. 

Default  value  of  0.85  will  be  used  to  select 
the  discharge  value  on  the  duration  curve 
used  to  calculate  dependable  energy.  Must 
be  zero  if  variable  DCAP  (PB.3)  is  positive. 

Streamflow  reliability  percentage  (decimal 
value  on  the)  used  to  select  discharge  value 
flow  duration  curve  for  use  in  calculation 
of  dependable  capacity. 

No  adjustment  will  be  made  to  the  energy 
benefit  determined  from  Table  1A  (page  12) 
of  text. 

Adjustment  factor  (decimal)  to  be  multiplied 
times  the  FERC  energy  values  (1978  dollars) 
from  Table  1A  of  text. 
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FIELD 

8 


VARIABLE  VALUE 


DESCRIPTION 


CRATIO 


SEBR 


0  No  adjustment  will  be  made  to  the  capacity 

benefit  obtained  from  Table  IB  (page  13)  of 
text. 

+  Adjustment  factor  (decimal)  to  be  multiplied 

times  the  FERC  capacity  value  (1978  dollars) 
from  Table  IB  (page  14)  of  text. 

0  Default  value  of  0.50  will  be  used  for  the 

adjustment  factor  to  be  multiplied  times  the 
firm  energy  benefit  vo  determine  the 
secondary  energy  benefit  rate  (see  EF  card). 

+  Adjustment  factor  (between  0  and  1.0)  to  be 

multiplied  times  time  firm  energy  benefit  to 
determine  the  secondary  energy  benefit  rate 
(see  EF  card) . 
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5.2  CB  CARDS  -  Capacity  Benefit  Cards 


The  CB  cards  are  provided  when  capacity  benefit  cards  are  required 
when  variable  IREG  (PB.1)  is  zero  (or  PB  card  is  omitted),  and  benefit 
calculations  are  desired.  The  capacity  benefit  is  provided  as  a  function 
of  the  annual  plant  factor  ( APF ) .  The  11  fields  (2  cards  required) 
indicate  the  benefit  corresponding  to  values  of  APF  from  0  to  1.0  in 
increments  of  0.10.  For  a  constant  vaue  of  capacity  benefits  (regardless 
of  plant  factor),  input  the  same  capacity  benefit  eleven  times. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

CB 

Card  identification  characters. 

1-11 

CB 

+ 

Capacity  benefit  corresponding  to 
appropraite  APF  in  units  of  $/kw-yr.  Field 

11  is  input  in  first  field  of  second  CB 
card. 

12 

CRATIO 

0 

No  adjustment  will  be  made  to  capacity 
benefits  provided  on  the  first  11  fields  of 
the  CB  cards. 

+ 

Adjustment  factor  (decimal)  to  be  multiplied 

times  the  capacity  benefits  on  Fields  1-11 
of  the  CB  card.  Field  12  is  input  in  the 
second  field  of  second  CB  card. 
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5.3  EB  CARDS  -  Energy  Benefits  Cards 


Energy  benefit  cards  are  required  when  variable  IREG  (PB.1)  is  zero 
(or  PB  card  - r  omitted),  and  benefit  calculations  are  desired.  The  cards 
indicate  the  energy  benefit  is  a  function  of  APF.  for  a  constant  value  of 
energy  benefit  (regardless  of  plant  factor),  input  the  same  energy  benefit 
eleven  times . 


FIELD  VARIABLE  VALUE 


DESCRIPTION 


0  ID  EB 

1-11  EB  + 


12  ERATIO  0 


Card  identification  characters. 

Energy  benefit  corresponding  to  appropriate 
APF  in  units  of  $/mwh.  Field  11  is  input  on 
the  first  of  the  second  EB  cards. 

No  adjustment  will  be  made  to  the  energy 
benefits  provied  on  the  first  eleven  fields 
of  the  EB  cards. 


+  Adjustment  factor  (decimal)  to  be  multiplied 

times  the  energy  benefits  specified  on  the 
first  eleven  fields  of  the  EB  card.  Field 
12  is  input  on  the  second  field  of  the  EB 
card. 
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MC 

ME 

5.4  MC  CARD { S )  -  Seasonal  Capacity  Benefits 


The  MC  card  specifies  seasonal  capacity  benefits.  The  use  of  the  MC 
card  overrides  values  input  for  variable  IREG  (PB.1)  or  the  CB  card. 


FIELD  VARIABLE  VALUE  DESCRIPTION 


0  ID  MC  Card  identification  characters. 

1-11  CMAF  +  or  0  Capacity  benefit  for  each  season  required. 

The  eleventh  value  is  input  in  the  first 
field  of  the  second  MC  card. 


5.5  ME  CARD ( S )  -  Seasonal  Energy  Benefits 


The  ME  card  specifies  seasonal  enegy  benefits.  The  use  of  ME  card 
overrides  input  values  for  variable  IREG  (PB.1)  or  for  the  EB  card. 


FIELD  VARIABLE  VALUE  DESCRIPTION 


0  ID  ME  Card  identification  characters. 

1-11  EMAF  +  or  0  Energy  benefit  for  each  season.  The 

eleventh  value  is  input  in  the  first  field 
of  the  second  ME  card. 
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5.6  EF  CARDS  -  Equivalent  Firm  Energy  Benefits 

The  EF  card  specifies  the  value  of  secondary  energy  in  terms  of  an 
equivalent  firm  energy  value.  The  equivalent  firm  energy  value  may  range 
between  0.0  and  1.0  and  is  related  to  the  percent  of  time  the  secondary 
energy  is  available.  As  the  availability  of  secondary  energy  approaches 
100  percent,  its  value  becomes  equal  to  the  firm  energy  value. 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

EF 

Card  identification  characters 

1-11 

SEFV 

+  or  0 

Equivalent  firm  energy  values  corresponding 

to  the  availability  of  the  secondary  energy 
in  terms  of  an  annual  percentage.  Input 
eleven  (11)  values  corresponding  to  the 
following  availability  percentages: 
(0,10,20,30,40,50,60,70,80,90.100).  The 
eleventh  value  is  input  in  the  first  field 
of  the  second  EF  card. 
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6  POWER  COST  INFORMATION 


Power  cost  cards  (Cl,  C2  and  CF  cards)  specify  cost  information  for 
analysis  desired. 


6.1  Cl  CARD  -  First  Card  -  Cost  of  Construction 


The  Cl  card  specifies  information  on  the  size  of  the  project  and 
related  components. 


FIELD  VARIABLE  VALUE 


DESCRIPTION 


0  ID 

1  HEIGHT 

(0) 


2  DIST 

(0) 


3  KS 

(0) 


Cl  Card  identification  characters. 

0  No  embankment  costs  are  calculated. 

+  Height  of  the  dam  in  feet  or  meters.  This 

value  is  used  in  the  computation  of 
embankment  costs. 

0  No  embankment  costs  are  calculated. 

+  Length  of  the  dam  crest  in  feet  or  meters. 

This  value  is  used  in  the  computation  of 
embankment  costs . 

Valley  shape  coefficient.  This  value 
is  used  to  determine  estimated  embankment 
costs .  Various  shapes  are  shown  Figure  1 
(page  38). 

0  No  embankment  costs  are  calculated. 

1  Trapezoidal  (Type  A) 

2  Triangular  (Type  B) 

3  Trapezoidal  river  section  with  a  trianglar 

over bank  (Type  C). 
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FIELD 

4 


5 


VARIABLE  VALUE 


DESCRIPTION 


ECAP 

(0) 


OPER 

(L) 


0  Existing  capacity  is  not  in  place. 


+  Existing  capacity  at  proposed  site.  The 

specified  or  calculated  value  of  CAPDES 
(PD. 4)  (installed  capacity)  is  assumed  to  se 
the  capacity  added  to  ECAP. 

L  Powerhouse  is  locally  operated.  Used  in 

estimating  OM&R  costs.  User  should  assume 
local  operation  unless  it  is  known  that 
project  is  remotely  operated. 

R  Powerhouse  is  operated  from  a  remote  site. 

Used  in  OM&R  cost  estimates. 


TURB  Three  character  code  for  setting 
(0)  extrapolation  flags  in  determining 

powerhouse  cost  and  for  indicating  the  type 
of  turbine  to  be  installed.  If  zero, 
extrapolation  of  powerhouse  costs  will  not 
occur.  The  program  will  select  the  types  of 
units . 


The  first  or  second  characters  may  be: 

C  Extrapolation  of  capacity  arrays  is  allowed 

in  determining  powerhouse  costs. 

H  Extrapolation  of  head  arrays  is  allowed  in 

determining  powerhouse  costs. 

NOTE:  If  neither  code  is  specified, 
extrapolation  is  not  allowed  and  a  total 
annual  cost  of  9999999999  is  assumed  for  the 
installation. 


The  third  character  (or  first  if  C  and  H  not 
used)  may  be: 

0  Program  selects  the  type  of  turbine. 

F  Francis  unit. 

K  Kaplan  unit  (usually  above  10  megawatts). 

S  Small  kaplan  unit  (less  than  10  megawatts). 

T  Tube  turbine 
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Cl 

FIELD  VARIABLE  VALUE  DESCRIPTION 

7  RESA  0  Reservoir  acquisition  and  clearing  costs 

(0)  are  not  calculated. 

+  Reservoir  area  in  acres  or  square 

kilometers . 

Used  to  determine  land  acquisition,  clearing 
and  preparation  costs . 

8  WYL  0  No  waterway  costs  are  calcualted. 

(0) 

+  Length  in  feet  or  meters  of  waterway  or 

diversion  canal  that  transports  water  to  the 
powerhouse . 

9  WYQ  0  No  waterway  costs  are  calculated. 

(0) 

+  Maximum  discharge  in  cfs  or  cms  associated 

with  the  waterway  specified  in  Cl. 8. 

10  IPROJ  Code  indicating  components  included  in  the 

(0)  cost  calcuations.  The  sum  of  the  following 

codes  deletes  the  specified  project 
components  from  the  cost  calculations.  When 
IPROJ  is  specified,  it  will  override  a 
request  for  cost  computations  on  Fields  1-9. 

0  All  components,  including  a  waterway,  are 

included  in  the  cost  calculations. 

1  Delete  waterway  costs.  (Fields  Cl. 8  and 
Cl. 9  do  not  need  to  be  specified.) 

2  Delete  embankment  costs.  (Fields  C1.1, 

Cl. 2  and  Cl. 3  do  not  need  to  be  specified.) 

4  Delete  spillway  costs. 

8  Delete  reservoir  clearing  costs. 

16  Delete  land  acquisition  costs. 

32  Delete  cost  of  the  inlet  and  outlets. 

64  Delete  costs  of  the  powerplant 

128  Delete  operation,  maintenance,  and 

replacement  costs. 
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FIGURE 


TYPE  A. 


TYPE  B. 


TYPE  C. 


:  VALLEY  SHAPE  CODES 


1 - 1 - -I 

-  xxxxxxxxxxxxxxxxxxxx 

lx  X 

h  x  x 

lx  X 

xxxxxxxxxxxx 


TRAPEZOIDAL  RIVER  CROSS  SECTION 


1 

1 

h 

1 


xxxxxxxxxxx 

X  X 

X  X 

X  X 
X  X 
X 


TRIANGULAR  RIVER  CROSS  SECTION 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


1  > 

h 

1 

1  h/2 


X  X 

X  X 

X  X 

xxxxxx 

-1/6— 


TRAPEZOIDAL  RIVER  CROSS  SECTION  WITH  TRIANGULAR  OVERBANK 
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6.2  C2  CARD  -  Second  Card  -  Cost  Construction 

The  C2  card  enables  users  to  specify  actual  rather  than  general  (Cl 
card)  cost  information.  Whenever  actual  costs  are  provided  they  will 
override  costs  that  are  automatically  calculated.  When  ECAP  (Cl. 4),  the 
existing  capacity,  is  greater  than  zero,  the  costs  in  the  fields  C2 
through  C9  represent  costs  to  be  added  to  the  existing  costs  for 
determination  of  total  costs.  All  costs  are  expressed  in  thousand  dollars 
(1000) . 

NOTE:  Fixed  cost  components  for  an  optimization  run  (variable  CAPDES 
(PD. 4)  is  greater  than  zero)  should  be  input  for  cost  items  that  are  not 
considered  a  function  of  installed  capacity. 


FIELD  VARIABLE  VALUE 


DESCRIPTION 


0  ID  C2  Card  identification  characters. 

1  JSTATE  +  Two  digit  code  from  Table  (page  27  of 

Exhibit  7)  indicating  the  state  the  project 
is  located.  Used  to  adjust  land 
acquisition,  construction  and  replacement 
costs. 


2  CIO 

(0) 


3  CPWH 

(0) 


0  Inlet  and  outlet  costs  are  automatically 

calculated. 

+  Cost  estimate  of  inlet  and  cutlet 

facilities . 

Adjustment  factor  to  be  appled  to  the 
computer  calculated  value.  A  value  of  -1.23 
would  increase  the  value  by  23  percent. 
NOTE:  When  determining  costs  for  adding 
capacity  to  an  existing  installation 

(variable  ECAP  (Cl. 4)  is  greater  than  0) 
this  value  represents  the  percentage  of  the 
original  cost  of  this  item  to  be  allocated 
to  the  hydropower  addition.  A  value  of 
-0.23  would  allocate  23  percent  of  the 
original  inlet  and  outlet  cost  to  the 

addition . 

0  Cost  of  powerhouse  is  automatically 

calculated. 

+  Cost  of  powerhouse. 
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FIELD 


4 


5 


6 


7 


8 


VARIABLE  VALUE 


DESCRIPTION 


CEMB 

(0) 


CSPW 

(0) 


CWWY 

(0) 


CRC 

(0) 


CLA 

(0) 


CMIS 

(0) 


Similar  to  explanation  for  cost  of  inlet  and 
outlet  (Field  C2.2). 

0  Cost  of  embankment  is  automatically 

calculated. 

+  Cost  of  embankment. 

Similar  to  explanation  for  cost  of  inlet  and 
outlet  (Field  C2.2). 

0  Cost  of  spillway  is  automatically 

calculated. 

+  Cost  of  spillway. 

Similar  to  explanation  for  cost  of  inlet  and 
outlet  (Field  C2.2). 

0  Cost  of  waterway  is  automatically 

calculated. 

+  Cost  of  waterway. 

Similar  to  explanation  for  cost  of  inlet  and 
outlet  (Field  C2.2). 

0  Cost  of  reservoir  clearing  is 

is  automatically  calculated. 

+  Cost  of  reservoir  clearing. 

Similar  to  explanation  for  cost  of  inlet  and 
outlet  (Field  C2. 2). 

0  Cost  of  land  acquisition  is 

automatically  calculated. 

+  Cost  of  land  acquisition. 

Similar  to  explanation  for  cost  of  inlet  and 
outlet  (Field  C2.2). 

0  No  miscellaneous  costs  are  added. 


+  Any  additional  fixed  costs  in  dollars  not 

provided  in  other  fields  of  the  C2  card. 
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FIELD  VARIABLE  VALUE 


DESCRIPTION 


10  COSTR 
(0) 


0  No  adjustment  factor  will  be  applied 

to  input  costs  estimates. 


+  Adjustment  factor  (decimal)  to  be  multiplied 

times  the  annual  cost  (Item  20  in  the  cost 
estimate  form). 
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6.3  CF  CARD  -  Cost  of  Financing  the  Project  Card 

The  CF  card  specifies  cost  information  necessary  to  calculate  the 
repayment  of  the  project  costs. 


FIELD  VARIABLE  VALUE 


DESCRIPTION 


0  ID 


CF  Card  identification  characters. 


1  CONT 

(0.25) 


0  A  contingency  factor  of  .25  (25  percent) 

will  be  used. 


+ 


2  RATE  0 

(0.06875) 


Contingency  factor  (decimal)  applied  to  the 
construction  costs.  A  value  of  0.25 
indicates  a  contingency  factor  of 
twenty-five  percent. 

An  annual  discount  rate  of  0.06b75  will 
be  used  in  the  computations. 


+  Annual  discount  rate  (decimal)  to  be  used  in 

the  computations. 


3  AMOP  0 

(100) 


Amortization  period  of  100  years  will  be 
used  in  the  computations. 


4  PTC 

(2) 


+  Length  in  years  of  the  amortization  period 

to  be  used  in  the  computations. 

+  Project  time  of  construction  in 

years.  Suggested  times  from  the  NPD  (1979) 
cost  estimating  manual  are  shown  below. 


Addition  of  a  power  plant  to  an  existing  dam 
-  2  years. 


Diversion  project  or  small  hydropower  and 
dam  project  (dam  height  less  than  100  feet) 

-  4  years .  Medium  hydropower  and  dam 
(project  dam  height  100  to  250  feet)  -  5 
years . 

Large  hydropower  and  dam  project  (dam  height 
greater  than  250  feet)  -  6  years. 

5  REPL  +  Replacement  cost  factor. 

(0.0125) 


EXHIBIT  7 


Page  45  of  54 


oc 


HYDUR  INPUT  DESCRIPTION 


7  OC  CARD  -  Optimization  Criterion  Card 


The  OC  card  specifies  the  information  necessary  to  define  the 
objective  function  used  to  optimize  the  installed  capacity.  Optimum 
capacity  selection  will  automatically  be  performed  if  CAPDES  (PD. 4)  is 
omitted.  If  the  OC  card  is  omitted  output  variable  26  (Table  2),  annual 
net  benefit,  will  be  maximized. 


FIELD  VARIABLE  VALUE 


DESCRIPTION 


0  ID 

1  IVAR1 

(26) 


2  IVAR2 

(0) 


3  ICOMB 

(0) 


OC  Card  identification  characters. 

0  Total  annual  net  benefits  ($/year)  will  be 

optimized,  (both  cost  and  benefit 
information  must  be  provided  on  the  PB,  CB, 
EB,  Cl,  C2,  CF  cards). 

+  Two  digit  numerical  code  indicating  what 

variable  in  Table  2  (page  16  of  40  of 
Exhibit  7)  is  to  be  maximized  or  minimized. 

0  Only  the  variable  specified  in  field  OC.1 

is  to  be  optimized. 

+  Two  digit  numerical  code  indicating  the 

second  variable  in  Table  2  to  be  combined 
VAR1 . 

This  field  specifies  what  type  of 
combination  of  VARl  and  VAR2is  to  be 
optimized. 

0  Optimize  VARl  only 

1  Optimize  the  sum  of  VAR2  and  VARl 

2  Optimize  the  remainder  when  VAR2  is 
subtracted  from  VARl. 

3  Optimize  product  of  VARl  and  VAR2. 

4  Optimize  ratio  of  VARl  and  VAR2. 
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oc 


VARIABLE  VALUE 


DESCRIPTION 


OPERND 

(MAX) 


This  field  indicates  whether  the 
optimization  function  is  to  be  maximized  or 
minimized. 

MAX  Select  the  installed  capacity  which 
maximizes  the  variable. 

MIN  Select  the  installed  capacity  which 
minimizes  the  variable. 
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8  OP  CARD  -  Optimization  Ordinate  Card 


The  OP  card  specifies  the  percent  of  time  exceeded  ordinates  used  in 
the  optimization  table.  The  default  percent  of  time  exceeded  ordinates  of 
1,5,10,20,40,60,80,90,95,  and  99  are  used  if  the  OP  is  omitted.  NOTE: 
This  card  is  only  operational  during  annual  power  analyses . 


FIELD 

VARIABLE 

VALUE 

DESCRIPTION 

0 

ID 

OP 

Card  identification  characters. 

1-10 

OP 

+ 

Percent  of  time  exceeded  ordinates  used  in 
the  optimization  table. 
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9  PS  CARD  -  Printout  Suppression  Card 


The  PS  card  enables  the  user  to  specify  desired  output.  The  complete 
output  will  be  printed  if  the  PS  card  is  omitted.  The  user  must  enter  the 
sum  of  values  to  indicate  which  of  the  described  items  are  to  be 
suppressed  in  the  printout. 


FIELD  VARIABLE  VALUE  DESCRIPTION 


0  ID 


PS  Card  identification  characters . 


1  PLOTS 

(0) 


2  CEFS 

(0) 


3  ECHOS 

(0) 


4  CALCS 

(0) 


0  No  plots  will  be  suppressed. 

1  The  streamf low-duration  plot  is  suppressed. 

2  The  total  capacity-duration  plot  is 
suppressed. 

4  The  existing  capacity-duration  plot  is 

suppressed. 

0  No  cost  estimating  forms  will  be 

suppressed. 

1  The  total  potential  capacity  cost  estimate 
form  will  be  suppressed. 

2  The  existing  capacity  cost  estimate  form 
will  be  suppressed. 

0  The  summary  of  the  input  variables  and 

program  options  printed  at  the  beginning  of 
the  job  will  not  be  suppressed. 

1  The  summary  of  input  variables  and  program 

options  printed  at  the  beginning  of  the  job 
will  be  suppressed. 

0  None  of  the  following  items  will  be 

suppressed. 

1  The  tabulated  ordinates  of  the  adjusted 

flow-duration  curve  will  be  suppressed. 
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PS 

FIELD  VARIABLE  VALUE  DESCRIPTION 

2  The  table  summarizing  the  optimization 

calculations  will  be  suppressed. 

4  The  table  summarizing  the  power  potential 

results  will  be  suppressed. 

5  IPRNTS  Seasonal  printout  indicator 

(0) 

0  Only  final  (annualized)  information  will  be 

printed. 

4095  Print  all  season's  output 

(or  greater) 

1  Print  1st  season's  output. 

2  Print  2nd  season's  output 

4  Print  3rd  season's  output. 

8  Print  4th  season's  output. 

16  Print  5th  season's  output. 

32  Print  6th  season's  output. 

64  Print  7th  season's  output. 

128  Print  8th  season's  output. 

256  Print  9th  season's  output. 

512  Print  10th  season's  output. 

1024  Print  11th  season's  output. 

2048  Print  12th  season's  output. 
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FIELD  VARIABLE  VALUE 


DESCRIPTION 


10  TRACE  0 

1 

2 

4 

8 

16 

32 

64 

128 

256 


Suppress  all  program  trace  output. 

Trace  power  calculation  for  exisitng  or 
incremental  capacity. 

Output  existing  poverarrays. 

Output  incremental  power  arrays. 

Output  optimization  power  arrays. 

Output  fixed  power  incremental  arrays. 

Output  fixed  power  arrays. 

Output  optimization  matrix. 

Trace  optimization  function. 

Trace  QGUESS/SELQ  calculations. 
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10  EJ  CARD  -  End  of  Job  Card  ( Required) 

This  card  must  be  provided  at  the  end  of  each  job.  Multiple  jobs  can 
be  executed  by  providing  subsequent  sets  of  T1  through  EJ  cards. 

FIELD  VARIABLE  VALUE  DESCRIPTION 

0  ID  EJ  Card  identification  characters. 
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